Prepared for submission to JHEP 



RBRC-863 
SHEP-1043 



One-loop operator matching in the static heavy and 
domain-wall light quark system with 0{a) 
improvement 



Tomomi Ishikawa,"'' Yasumichi Aoki, Jonathan M. Flynn,'^ Taku Izubuchi^'^ and 
Oleg Loktik^ 

Physics Department, University of Connecticut, 
Starrs, CT 06269-30^6, USA 

^ RIKEN BNL Research Center, Brookhaven National Laboratory, 

Upton, NY 11973, USA 
'^School of Physics and Astronomy, University of Southampton, 

Highfield, Southampton SO 17 IB J, UK 
^'■Physics Department, Brookhaven National Laboratory, 

Upton, NY 11973, USA 
'^Physics Department, Columbia University, 

New York, NY 10027, USA 

E-mail: tomomi@phys.uconn.edu, yaoki@quark.phy.bnl.gov, 
j . m . f lynn@soton . ac . uk, izubuchi@quark . phy . bnl . gov, 
oleg . loktik@gmail . com 

Abstract: We discuss perturbative 0{g^a) matching with static heavy quarks and 
domain-wall light quarks for lattice operators relevant to i?-meson decays and B^-B^ 
mixing. The chiral symmetry of the light domain-wall quarks does not prohibit opera- 
tor mixing at 0(a) for these operators. The 0{a) corrections to physical quantities are 
non-negligible and must be included to obtain high-precision simulation results for CKM 
physics. We provide results using plaquette, Symanzik, Iwasaki and DBW2 gluon actions 
and applying APE, HYPl and HYP2 link-smearing for the static quark action. 

Keywords: Lattice QCD, Heavy Quark Physics 



^Present address: Kobayashi-Maskawa Institute for the Origin of Particles and the Universe (KMI), 
Nagoya University, Nagoya 464-8602, Japan. 



Contents 



1 Introduction 1 

2 General discussion of matching in HQET 2 

2.1 Definition of the heavy quark field and state normalization 2 

2.2 Operators relevant for B meson decays and — B mixing 3 

2.3 Matching procedure 4 

2.4 HQET operators and symmetries 5 

2.4.1 Symmetries in the static heavy and domain- wall light quark system 5 

2.4.2 Quark bilinear operators 6 

2.4.3 AB = 2 four-quark operator 7 

2.5 Continuum matching: CQCD i — > CHQET 9 

2.6 Static effective theory matching: CHQET i — > LHQET 12 

3 One-loop perturbative lattice to continuum matching in HQET 14 

3.1 Quark bilinear operators 14 

3.2 AB = 2 four-quark operator 17 

4 Lattice perturbation theory 18 

4.1 Lattice action 18 

4.1.1 Standard static action with link smearing 19 

4.1.2 Domain-wall fermion action 19 

4.1.3 Gluon action 20 

4.2 Mean field improvement 21 

4.3 Alternative form for the 0{pa) terms 24 

4.4 Power divergence from 0{pa) operator 25 

5 Numerical value of the matching factor: an example 27 

5.1 Continuum matching factor 27 

5.2 Lattice to continuum matching factor and 0(a) improvement coefficient in 
HQET 28 

5.3 0(a) corrections for physical quantities 30 

6 Conclusion 31 

A Lattice Feynman rules 31 

B Calculation of lattice perturbation 34 

C Relations of the Feynman rules and the one-loop amplitudes betw^een 

different domain- wall fermion actions 38 



D Power divergence structure in equations of motion 40 



E Operator mixing in the static heavy and Wilson Ught quark system 



42 



F Tables of numerical values of domain-wall light quark part 



45 



G Tables of numerical values of integrals 



48 



1 Introduction 

The Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix elements [1, 2] play a key 
role in elementary particle physics. Constraints on Vts and Vtd can be obtained from 
mixing, where the 5C/(3)-breaking ratio ^ = /b^ \JBbJ Jb^ ^/Bb^ [3] plays an 
important role. Lattice QCD simulations for S-meson physics are challenging, however, 
because of the difference in energy scales between the light u and d quarks and the heavy 
b quark. One approach to this issue is to use Heavy Quark Effective Theory (HQET) [4]. 
Lattice calculations with HQET have several known difficulties (and solutions), (i) The 
static propagator is noisy because the static self-energy contains a 1/a power divergence. 
The ALPHA collaboration introduced a modified static action which improves the signal 
to noise ratio (S/N) [5]. The modification is to replace the gluon link variable in the 
static action with a smeared one obtained by 3-step hypercubic blocking. Using this action 
the power-divergent contributions in the static self-energy are significantly reduced, (ii) 
Non-perturbative matching to the continuum is needed. If we include 0{\/mh) corrections 
in HQET, the continuum limit cannot be reached using perturbative matching because 
of a power divergence [6]. Non-perturbative matching schemes include the Schrodinger 
functional with step scaling or RI/MOM, but implementing them in practice is not easy 



The static approximation (lowest order of HQET) is theoretically simple and its im- 
plementation in lattice calculations is relatively easy. Because of this it is often used in 
simulations as a first step. This leads to errors of ©(Aqcd/""^;)) ~ 10%. Although 0(l/mf,) 
corrections should be included for more precise calculations, the static approximation works 
well in the determination of the SU (3)-breaking ratio ^, in which the theoretical uncer- 
tainty is suppressed by (m^ — md) / ^.q^cd and is estimated to be about 2%. In addition, 
perturbative matching is more justified for the static approximation. Hence the aim of this 
paper is to calculate the one-loop perturbative matching factor for the static approximation 
to enable the determination of ^. 

Discretization errors in lattice simulations can limit the precision with which physical 
quantities are determined. Therefore the 0{a) improvement program on the lattice is 
important. The 0{a) improvement of heavy-light currents with clover Wilson light quarks 
was investigated using one-loop perturbation theory in non-relativistic QCD [7] and the 
static approximation [8]. The 0{a) effects give large corrections to B meson quantities. 
In this work we consider the case where the light quarks are simulated using a lattice 
action with good chiral symmetry. In the light-light quark system, the chiral symmetry 



for HQET. 
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guarantees the absence of 0{a) errors in the operator. For the static heavy-hght quark 
system, however, there are 0(a) effects even if we use chiral fermions for the hght quarks. 
This result was aheady found for the clover Wilson light quark with Wilson parameter 
r = (which is chirally symmetric, but has doublers) [8]. In this study we use the domain- 
wall (DW) fermion action [9-11] to realize chiral symmetry for the light quarks. This paper 
is an extension of several previous works: the perturbative matching at 0{g'^) in the static 
heavy and domain-wall light quark system [12] and its version with a link-smeared static 
action [13, 14]. The main results are the matching between continuum HQET and lattice 
HQET and 0(a) improvement of the heavy-light quark bilinear operator (Eq. (4.49)) and of 
the AS = 2 four-quark operator (Eq. (4.50)). The calculations are performed using mean- 
field (MF) improved one-loop lattice perturbation theory and include the link smearing of 
the static action. 

The paper is organized as follows. In Section 2 we give our field notation and discuss 
the matching procedure. We give symmetry-based restrictions on the forms of the HQET 
operators; these apply for non-perturbative treatments as well. In Section 3 we give the 
form of the one-loop matching expressions while in Section 4 we provide more discussion 
of the lattice perturbation theory calculations. Expressions for the continuum to HQET 
to lattice HQET matching with mean-field improvement are given in Eqs. (4.49) to (4.51). 
In Section 5 we compute numerical values of the one-loop matching factors and estimate 
the 0(a) effect for physical quantities including the SU{2>) breaking ratio ^. Concluding 
remarks are given in Section 6. The lattice Feynman rules are listed in Appendix A and 
details of the lattice perturbation theory calculation are shown in Appendix B. Numerical 
values of several integrals are listed in tables in Appendices F and G. 

2 General discussion of matching in HQET 

2.1 Definition of the heavy quark field and state normalization 

We regard the h quark as heavy and give it an on-shell velocity v = (1,0,0,0); then the 
on-shell momentum is p = rrn,v = (m;,, 0,0,0), where mi, is the b quark mass. We build a 
heavy quark field h as a sum of quark, and anti-quark, 

h{x) = h+{x) + h-{x), h{x) = h+{x) + h-{x) = h\{x)-h}_{x), (2.1) 

where 

h±{x) = e»'^-^i^6(a;) = e^"'''*i^^6(x). (2.2) 

In the static limit, /i+ and h- decouple. The definition of the creation and annihilation 
operators for them is as follows: 

h}j^ : creates the outgoing heavy quark, 

cannot annihilate the incoming anti-heavy quark. 

/i+ : annihilates the incoming heavy quark, 

cannot create the outgoing anti-heavy quark. 
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h- : creates the outgoing anti-heavy quark, 

cannot annihilate the incoming heavy quark. 

hl_ : annihilates the incoming anti-heavy quark, 
cannot create the outgoing heavy quark. 

The HQET state normalization differs from that normally used in QCD. The relation 
of the states is 

|5>QCD = V^[|5>HQET + 0(1M)], (2.3) 

where is the mass of the state |i?)QCD- 

We use the PDG notation in which the quark content of the B meson is B = (bq) and 
B = (bq), where q denotes light d or s quarks [15]. 

2.2 Operators relevant for B meson decays and B^ — B^ mixing 

The QCD operators considered in this paper are the heavy-light quark bilinear 

^QCD ^ y^^^ ^2.4) 

where T = {1, 75, 7^*75, (^iJ.u] with a^u = ^[ln,lu], and the AS = 2 four-quark operator 

oT'' = H<i\H^]^ (2-5) 

where 7^ = 7^Pr = 7/,(l + 75)72 and 7^ = 7^Pl = 7^,(1 - 75)/2. Quark colors are 
contracted within the square brackets. The static effective theory operators needed for 
matching to the QCD operators are 

^HQET ^ -j^^^^ ^2.6) 

OT""^ = [hlf:q][hll:q], (2.7) 

,H( 
S 



OjQ^^ = [hPLq][hPLq]. (2. 



Omitting parts of these operators which do not contribute in the matrix element for B^ — B^ 
mixing, Eqs. (2.7) and (2.8) become 

= [h+jf:q][h.jf:q] + [h.^l:q][h+^f:q] = 2[7I+7^g][7^_7^g], (2.9) 
O^'^''^ = {hPLqjlhPLqj 

= [h+PLq][h-PLq] + [h-PLq][h+PLq] = 2[h+PLq][h-PLq], (2.10) 

where we define the doubled square bracket by = [/i4_X][/i_y] -|- [/i_X][/i+y]. 
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2.3 Matching procedure 

We adopt a two step matching procedure: 



continuum QCD 




' - 

continuum HQET 






(CQCD) 


fi = nib 


(CHQET) 










RG-evoIution 






continuum HQET 




lattice HQET 




(CHQET) 


fi = a ^ 


(LHQET) 



in which we first perform the matching between continuum QCD (CQCD) and continuum 
HQET (CHQET) and subsequently match CHQET to lattice HQET (LHQET). Some 
comments on this matching: 

(1) The continuum QCD (CQCD) operators are renormalized in MS(NDR) at scale /Xf, 
which is usually chosen to be the b quark mass mb- Fierz transformations in arbitrary 
dimensions are specified in the NDR scheme introduced by Buras and Weisz [16]. The 
introduction of evanescent operators gives vanishing finite terms at one-loop but is 
needed to obtain the correct anomalous dimensions at two-loop. 

(2) The CHQET operators are also renormalized in MS (NDR) at some scale //. Matching 
between the continuum theories is performed in perturbation theory by calculating 
and comparing matrix elements of the operators between an initial state |i) and 
final state |f) for each theory. The calculation has been done for quark bilinears at 
one-loop [4] and two-loop [17] levels, and for the AB = 2 four-quark operator at 
one-loop [18]. 

(3) The continuum matching between QCD and HQET is done at scale /_i = nib to avoid 
a large logarithm of ^/nib- We use renormalization group (RG) running in CHQET 
to move to a lower scale at which the HQET matching between continuum and lattice 
is done. We employ the two-loop anomalous dimension calculations in Refs. [19, 20] 
for the bilinear and in Refs. [21-23] for the four-quark operator. 

(4) Matching between CHQET and LHQET is performed at scale /x = a~^, where a 
denotes the lattice spacing. The calculation is performed in one-loop perturbation 
theory taking into account 0{a) discretization errors on the lattice. For this we 
introduce external momenta, e.g., 

(f|0|i) = {h{p')\0\q{jp)) for bilinear operators, 

(f|0|i) = {h{p'2) q{p2)\0\h{p'i) , q{pi)) for four-quark operators, (2.11) 

where O denotes the bilinear and four-quark operators. On-shell improvement is 
used, in which we impose the equations of motion on the external quarks: 

1)0/1 = 0, (^ + mg)g = 0, (2.12) 
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where niq denotes the light quark mass. The transition amphtudes are expanded 
around zero in the external momenta up to first order, giving rise to 0{pa) effects. 
The 0{p'a) effects always vanish due to the on-shell condition p' = 0. The amplitudes 
are also expanded to first order in the light quark mass to reveal 0{ma) effects. 

2.4 HQET operators and symmetries 

In this section we introduce symmetry transformations that appear to be powerful tool to 
study the mixing structure of the HQET operators under consideration. The arguments in 
this section are used in Sees. 2.5 and 2.6. Similar discussions were presented in Refs. [24, 25]. 

2.4.1 Symmetries in the static heavy and domain- wall light quark system 

We use following symmetries to restrict the operator mixing. 

Chiral symmetry 

We assume an SU{Np) SU{Nf) chiral symmetry for the light quark sector (Np denotes 
number of light quark flavors) : 

QR — > Unqn, QL — > Ulql, (2.13) 

where qr = Prq, ql = PlQ, and Ur and Up represent the SU{Np) chiral transformation 
matrices. To account for explicit symmetry-breaking by the light quark masses we assume 
that the light quark mass matrix M = diag(mi, • • • ,m]\fp) transforms as follows: 

M UpMUl, UrM^uI, (2.14) 

This restricts the 0{ma) operator mixing. We use the domain- wall fermion formulation 
for light quarks on the lattice in this paper and treat these quarks as exactly chiral. In 
practice the domain-wall quark action breaks chiral symmetry by an amount proportional 
to powers of the residual quark mass which is of O(10~^) and hence small. 

Heavy quark spin symmetry (HQS) 

The static heavy quark obeys the heavy quark symmetry [26, 27] 



h± — Vih±, h± — h±vj, with Vj = exp 



(2.15) 



where A:} G {1,2,3} and (f)i is a rotation parameter. This symmetry comes from the 
spin-independence of the static quark interaction. 

Spatial rotational symmetry 

In the static limit the 4-dimensional Lorentz symmetry (hypercubic symmetry on the lat- 
tice) breaks to 3-dimensional rotational symmetry. For this reason we treat the temporal 
and spatial directions separately. 
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Discrete symmetries 

The static and domain-wall actions respect the discrete symmetries: parity V, time reversal 
T and charge conjugation C. The transformations are given by (we follow the definitions 
in the pedagogical reference [28].): 

V : q{x, t) — > -iQq{-x, t), q{x, t) — > q{-x, t)-fo, 

h±{x,t) — > ±h±{-x,t), h±{x,t) — > ±h±{-x,t), (2.16) 

T : q{x, t) — > 7o75q{x, -t), q{x, t) — > q{x, -t)757o, 

h±{x,t) — > ±75/izp(x, -t), h±{x,t) — > ±7izp(x, -t)75, (2.17) 



C : q{x,t) — > Cq{x,t)'^ , q{x,t) — > —q{x,t)^C ^, 

h±{x,t) — >Ch^{x,tf, h±{x,t) — > -h^{x,t)^C~^, (2.18) 

C = 7072, Cjf,C~^ = -7J. 
Their combinations are also useful: 

V-T: q{x,t) ^-f5q{-x,-t), h±{x,t) ^h^{-x, -t)-/^, (2.19) 

T-C: q{x, t) 7075Cg(x, -tf , h±{x, t) Th±{x, -tfC'^-f^, (2.20) 
2.4.2 Quark bilinear operators 

We introduce higher dimensional operators to improve J±r at 0{pa) and 0{ma): 

J±rD = h±r{'y • B)g, J±rM = mJi^Tq = nigj^^, (2.21) 

where we use spatial rotational symmetry and equations of motion (following the on-shell 
improvement program) to reduce the number of operators. In the following we see how the 
HQET operators transform under the symmetries listed above in Sec. 2.4.1. 

- V 

The bilinear operators transform under parity as 

J±r ^ GJ±r, J±vD — > GJ±rD, J±vm ^ GJ±ym, (2.22) 
where G is defined by 7or7o = GT and takes the values +1 or —1. 

- HQS 

Under the HQS transformation the bilinear operators form multiplets classified by G. 
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- Chiral symmetry 

Chiral transformations interchange multiplets classified by G = ±1. The rotational 
direction in this transformation of J±td and J±tm is opposite to J±r- 

- V-T 

Under the V ■ T transformation, 

VT V-T V-T 

J±r > J^r, J±rD > J^td-, J±vm > J^vm- (2.23) 

- T-C 

Under the T • C transformation, 

J±r ^ -H^K ( J±r)^ , J±rD ^ -H^K (J±rz))^ , J±rM ^ -H^K ( J±rAf)^ , 

(2.24) 

where and K are defined by 75r75 = H^T and C^^TC = KT* respectively. 
2.4.3 AB = 2 four-quark operator 

We introduce dimension seven operators to improve the AB = 2 operators. Using the 
equations of motion, spatial rotational symmetry and chiral symmetry, the independent 
operators can be reduced to 

Q±ND = 2[7^±7;?(7 • ^)q][hTlf^q], Q'±nd = 'i[h±PRh ■ ^)q][h^PLq], (2.25) 
Q±NAI = 2mg [h±^l!q] [h^^llq] , Q'±nm = ^±PRq] [h^PLq] , (2.26) 

where Q±nd and Q'^-j^d are 0{pa) operators, Q±nm and Q'j_j^j^.j are 0{ma) operators. 
The factor 2 in Eqs. (2.25) and (2.26) is chosen for convenience. 

-V-T 

Under the V ■ T transformation, 

Ol ^ Ol, (2.27) 
-p. 7" ^ V-T f V-T I V-T I 

Q±ND > Qtnd, Q±nd ^ Qtnd^ Q±nm > Qtnai, Q±nm ^ QtNM- 

Thus the 0{pa) and 0{ma) operators are written as 

QnD = Q+ND + Q-ND, Q'nD = Q+ND + Q'~NDj (2.28) 

QnM = Q+NM + Q-NM, Q'nM — Q'+NM + Q'-NM- 

- T-C 

Under the T - C transformation, 

OL^(Oi)^ Qnd^{Qnd)^, Qnm^{Qnm)\ (2.29) 
Os (Os)^ , Q'nd [Q'nd]'^ I Q'nm [Q'nm)'^ - 
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HQS 



We define the following bilinear operators: 



(2.30) 



for convenience. Under an HQS transformation along the i-th direction (i-HQS), these 
operators transform as 

^ i~HQS J J i-HQS J ^ i-HQS ^ 
-'±0 > —J±i-, ^±1 > ^±0, J±j(^i) > ^ijkJ±k- 

i-HQS i-HQS i-HQS 

^±0 > J^±i, J^±i > --'■^±0, J^±j{:^i) > ^ijkJ^±k, 

where X ^^^> Y means Vi{X) = 1 + 4)iY + 0{(f)f). Next, we consider the four-quark 
operators: 



(2.31) 



where the free index a is not summed. Under the i-HQS transformation, these operators 
transform as 



Lq. = 




— J +aJ —a ~\~ J—aJ+a: 








Ra = 









i-HQS 


Li 


i-HQS 
> Lo, 




j-HQS 

> Lk{^i,j), 


No- 




Ni - 


™> iVo, 




^ Nk{^i,j) 


Ro 


i-HQS „ 
> -Kj, 


Ri 


i-HQS „ 
> -Ko, 




i-HQS „ 
^ ^K^iJ) 



(2.32) 



Using this, we can define the ^-HQS transformation as 
h± — V/i±, h± 



E-HQS,,, , E-_HQS^^^^t, withV = nV. 



(2.33) 



1=1 



such that 



Lr 



E-HQS 



Lj, Li 



E-HQS 



Rq 



E -HQS^ 



, Ri 



E-HQS 



Lo + '^ Lk, 

kj^i 

> -No + Y,^k, 

k^i 

Ro + ^ -Rfc. 



(2.34) 



Under the ^-HQS transformation. 



Vl 



Lo + E^Li 

-3Lo + E^L^ 



Vn 



-No + ZiNi 
3No + EiNi 



+ Ri 

-3Ro + Yli Ri 



(2.35) 
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behave as eigenvectors: 



E-HQS 

yL,N,R > 



3 
-1 



^L,N,R- 



(2.36) 



Applying the transformation behavior (2.36) and imposing the chiral symmetry, we find 
0(1), 0{pa) and 0{ma) eigenvectors which have the same symmetry transformation 
property as O^: 



Ol 
Ol + 40s 



Ond 



Onm 

o- 



'NM 



where we define: 

Ond = Qnd + 2Q'nd = 2[hlj!h ' ^)Q\[hll:q] + 4[7^Pk(7 • ^)q][hPL 



^ND — 

Onm = Qnm + 2Q^m = {2[h-fl!q][h-/f^q] + A[hPnq][hPLq]) 

2[hj;!q][h7f^q]-^hPRq][hPLq]) 



o- 



NM 



Qnm — 2Q/ 



(2.37) 



(2.38) 
(2.39) 
(2.40) 
(2.41) 



When the chiral symmetry is not imposed, more operators arise besides those in Eq. (2.37). 
This is the case both for the Wilson-type quarks, also for the domain-wall fermion with 
finite extent of the fifth dimension. We will discuss this issue in Appendix E. 

2.5 Continuum matching: CQCD i — > CHQET 

The operator matching between QCD and HQET in the continuum is expressed by 



O 



^Q^°(/.,) = Cr(^6,/^)Jr'''^''^(^) + O(AqcdM), (2.42) 

CHQET (^^^ + Z2(m6, ^^)0™^(^) + 0(AQCD/m,). (2.43) 



CQCD 



{^lb) = Zi{i^ib,fi)o^ 



We employ one-loop perturbative matching. For the bilinear operators, we quote Ref. [4]: 



Cr(^.,/i) = l + (|^)' 



Cf 



In 



(2.44) 



with 



Ar = CF 


(^Dqcd — -Dhqet + 


-E^QCD - -EhQET \ 

2 


(2.45) 


Dqcd = - 


^^ + ^H' HH' 
2 4 


- 1, -E'qcd = -4, 


(2.46) 




-Chqet 


= 1; -Ehqet = 0, 


(2.47) 



where Cf = {N'^ — l)/(2-^c) (-^c = 3 is the number of colors) is a second Casimir, and we 
introduced H defined by HT = ^^T^^ and H' which is the derivative of H with respect to 
d in d dimensions. For the four-quark operator, we quote Ref. [18, 23]: 



-61n( 1^ ) +Bl 



^2(Atfe,At) 



(47r) 



(2.48) 
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where 



(2.49) 



To avoid large logarithms of ^i/ ^ib in Eqs. (2.44) and (2.48) we match at scale fi = = 
rrib and use renormalization group running in the effective theory to reach a lower scale 
fjL < mb'. 



1^ 



_d_ 



1 



Cr(/xb,/i) = -Cr{fib,t^h] 



2 ^ 



Zl{fJ.b,l-L) ^2(^6,^) 



711 712 
721 722 



(2.50) 
(2.51) 



where 7's denote the anomalous dimension of each operator. Here HQS (2.36) ensures that 

712 = 0, 722 = 711 + 4721- (2.52) 
Solving the RG-equations, we obtain the RG-evolution of the renormalization factors: 



Cr{mb,iJ,) = Cr{mb,mb)Urimb,n), 
Zi{mb,fi) Z2{mb,fJ.) = Zi{mb,mb) Z2{mb,mb) 



U'^^\mb,i2) U'l'^\mb,lj) 
Uf^\mb,n) U^^'^\mb, fi) 



(2.53) 



(2.54) 



where the factors U{mb, fJ^Ys represent the RG-evolution on the HQET side, and the sym- 
metry constraints (2.52) give 

U^^^\7nb,fi) = 0, uP{mb,fi) = ui^^\7nb,fi) + mP{mb,fi). (2.55) 

The one-loop matching requires a two-loop calculation of the anomalous dimension 

2 



where Os = /{A-k), and of the beta-function 



/3(«s) 



a 



with 



/3o 



11 



34 



Pi + 0{al) 



Ca - -TpNf, /3i = -Ci -A{-Ca + Cf TpNf, 



(2.56) 



(2.57) 



(2.58) 



where Ca = -^c and Tp = 1/2 are color factors for the fundamental representation of 
SU{Nc) and Nf denotes number of flavors of quarks with masses below fi. The two-loop 
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anomalous dimension was calculated, in Refs. [19, 20] for the bilinear and in Refs. [21-23] 
for the four-quark operator. Here we quote the results: 



where 



4vr J [ UsifJ.) 



dr 



+ 0(a2), 







di 


'OsifJ-bY 


d2\ 













+ 0{as), 



as{^J'b) 



d2 



+ 0{as), 



(0) 



111 
2/3o 



(0) 
722 



(2.59) 
(2.60) 
(2.61) 
(2.62) 

(2.63) 



with 



7?' 



(0) 

Til 



111 



Jr 




, Ji 




-3Ci?, 



-Cf I^Ca - \cf - ^TpNf - A{Ca - 4Ci.)C(2)| 



(0) 
722 



-6Cf + 4 



iVc + 1 



1 

6 

+- 



Cf {UQCf - 190Ca + 56TFNf - 48CaC(2)} 
A^.-l 



(2.64) 

(2.65) 
(2.66) 

(2.67) 
(2.68) 



2iV. 



{-240Cf + 104Ca + (24 - 16TF)Nf + 48(8Cf - Ca)C{2)} 



where C(2) = it^/6. 

We assume the lattice cutoff is higher than the charm quark mass rric- Since dynamical 
lattice QCD simulations are usually performed with up, down and strange sea quarks, not 
including charm, we employ a two step running [29]: 



U 



(11) 

L 

r(21) 



Ur{mb,a ^) 

(22), 



(2.69) 



(12), 



(22), 



Nf 



u 



(11), 



rric, a 



(m„a-i) 



C/t'')(m„a-i) i7i'')(m„a-i) 



Nf 



, (2.70) 
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in which numerical values of Eq. (2.64) are given by 

J^f = 0.910, ^ = 1.864, (2.71) 
J^f=^ = 0.755, jf^'"^ = 1.698. (2.72) 

2.6 Static effective theory matching: CHQET i — > LHQET 

In the matching between continuum and the lattice HQET we include 0{pa) and 0{ma) 
contributions, where p and m represent the light quark external momenta and light quark 
mass respectively. There is no 0{a) correction for the light-light quark system with exact 
chiral symmetry, but this correction is not prohibited for the static heavy-light case. 

Quark bihnear 

To obtain the matching factors, we calculate the on-shell transition amplitude of the op- 
erator J±Y with initial state |i) and final state |f) with external momenta, as described 
in Sec. 2.3. Taking into account the 0{a) error, we define a column vector of the HQET 
operators 

1 T 



J±T aJ±rD aJ±rM 



(2.73) 



Using them, the transition amplitudes in the CHQET through the first order in 0{pa, ma) 
are given by 

(f|J±r|i)cHQET = U(^) GA'-P''^ ((f|J±r|i))cHQET 

L J cont 

= ^G,cont((f|J±r|i))cHQET, (2.74) 

where ((•••)) denotes a tree- level amplitude and we use the abbreviation [Ax Bx ■ ■ ■] = 
[A B • • • ]x. In this expression, are all real and independent of ±T. 

That fact and the existence of G in Eq. (2.74) are read from the transformation properties 
in Sec. 2.4.2. The transition amplitudes for the LHQET are generically written as 



(f|c7±r|i)LHQET 



((f|J±r|i))LHQET 



QJ[{pa;l) ^(pa.pa) 
^^(■ma,l) j^{'ma,pa) j^{ma,ma) 

J latt 

^G,latt((f|:4r|i»LHQET. (2.75) 



Here we have introduced the operator mixing matrix elements A^^^'^^ and A^^""'^^ which 
represent 1/a power-divergent mixing with lower dimensional operators. This possibility 
needs to be taken into account because the reduced symmetry in lattice theories may not 
be sufficient to prohibit power divergences [30]. We will discuss this more in Sec. 4.4. 
Combining Eqs. (2.74) and (2.75), the matching between the continuum and the lattice 
HQET can be written as 

(f| J±r|i)cHQET = ^G,contdiag[cS(J±r),5(aJ±rD),'5(aJ±rAf)] (-^^latt)"^ (f|.7±r|i)LHQET 

= .^r(f|</±r|i)LHQET + ^rDG(f|aJ±rD|i)LHQET + ^r7\/G'(f|aJ±rA/|i)LHQET, 

(2.76) 
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where S{X) = ((f|X|i))cHQET((f|-^|i))LHQET' '^^lich absorbs the difference in the tree-level 
amplitude between CHQET and LHQET. 

Ag and Aq^ can be expanded in perturbation theory both in the continuum and on 
the lattice as 



Ag = I+(^) CpAG + Oig'), Aa' = I-(^) CMg + 0(/). 



4-kJ 

At one-loop it is sufficient to calculate 



\4ttJ 



(2.77) 



(f|'^±r|i)cHQET 



Cf (41 - A^;;;^ ) } 5(J±r)(f| J±r|i)LHQET 



+ 



+ 



(- 

V47r 

-V 



CpG ( A^^2t 



CfG (4^,^ - A\t^r' ) 5(aJ±rM)(f|aJ±rM|i>LHQET + 0(/). (2.78) 



Air^)cS(aJ±rD)(f|aJ±rD|i) LHQET 



I (l,ma) 



Four-quark operator 

From Eq. (2.36) there are two distinct multiplets associated with the S-HQS transforma- 
tion. Defining column vectors of the HQET operators: 

1 T 



02 



Ol + 4O5 aO^ aOjjj^ 
the transition amplitudes in the CHQET can be written as 

(f|OL|i)cHQET 



(2.79) 
(2.80) 



^ J cont 
= ^1, cont((f|C'l|i))cHQET, 



|Cl|i))cHQET 



(f|Oz, +405|i)cHQET = 

= ^2, cont 

The transition amplitudes for the LHQET are generically 



io2iiy 



cont 
C2|i))cHQET- 



CHQET 



(2.81) 
(2.82) 



(f|0 

o I i) LHQET 



B, 



(1,1) 



„(l,pa) „{l,ma) 



<o{pa,l) y2{pa,pa) y,(pa,ma) 

l^a Oa Oa 

j3(m.a,l) yj{m,a,pa) yj{m,a,ma) 

l^a l^a 



((f|0„|i))LHQET 



latt 



B, 



a, latt\ 



LHQET, 



(2.83) 



where a takes on the values 1 or 2. Using these relations, the matching between the 
continuum and the lattice HQET can be written as 

(f|OL|i)cHQET (2.84) 
= i3Lcontdiag[5(OL),cS(aOArD),5(aOArAf)] (6i,iatt)~^ (f|Cl|i)LHQET 
= -^L(f|OL|i)LHQET + -^AfD (f|«OAf_D |i) LHQET + -^ArA/(f |«07VA/|i)LHQET, 

(f|(OL + 405)|i)cHQET (2.85) 
= B2, contdiag[5(0L + 40s),5(aO^),cS(aOjvA7)] (^2,iatt)"^ (f|02| 1/ LHQET 
= ZL+is{^\{OL + 40s)|i)LHQET + -^ivD^^ I'^^IvdWlhqet + Z^^{i\aOj^\\)u:{QiErr- 
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In perturbation theory both in the continuum and on the lattice, Ba and B^ ^ are expressed 
as 

2 . , , / n \ 2 



and such that expression for the matching between CHQET and LHQET is: 



(f| 



Ol 
Os 



/CHQET 





' Bi O" 


(1) 


' Bi O" 






cont 




latt / J 



Ol 
Os 





' Bi 




(pa) 




" Bi 




(l,pa)\ 












( 


Bi 








Bi 








4 


4 _ 


cont 




4 


4 _ 


latt / 



xdiag[5(OL),5(05)](f| 



X diag [S {aOND ) , 5 {aOjj^)] (f | 



/LHQET 



qOnd 
aO^ 



/LHQET 









(ma) 






(l,ma)\ 




' Bi 





" Si 





( 


Bi 


B2 












4 


4 _ 


cont 


4 


4 _ 


latt / 



xdiag[5(aOArAf),5(aO- 



cOnm 

aa 



'NM 



(2.86) 



|i)LHQET + 0(/).(2.87) 



3 One-loop perturbative lattice to continuum matching in HQET 

In the continuum perturbative calculations we use the MS renormalization scheme. Feyn- 
man gauge is employed and IR divergences are regulated using a gluon mass A in the 
continuum and lattice calculations. In this section we only collect the results. The actual 
calculations are shown in Sec. 4 and Appendices A and B. For simplicity, the lattice spacing 
a is set to be 1 in the calculations. 



3.1 Quark bilinear operators 

We consider one-loop perturbative matching of HQET operators between the continuum 
and the lattice. The one-loop transition amplitudes (2.74) and (2.75) for the HQET can 
be written as 



(f|J±r|i>HQET = |l+ (£)'cM(^)|((f|J±r|i))HQET 



(3.1) 



+ (|^)'cpG'i(P'^)((f|aJ±rD|i»HQET + (£)' C^Gi(™'^)((f|aJ±rA/|i))HQET. 
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(- ^ + ^HQET^ 




-^CHQET ~ 

l(ma) _ 
-^CHQET ~ 






3aX ' 
4tt 

aA' 





In this expression the continuum coefficients -^^^Qgrp, ^hqet and ^I^'^et are 

2 In ^ + Ehqet) + ^ (hi ^ + ^) , (3.2) 

(3.3) 
(3.4) 

where -Dhqet and -Ehqet are shown in Eq. (2.47) and 

F = i (3.5) 

The lattice coefficients are 

4'i1qet= (-ln(a^A2)+d«)+^(-21n(a^A2)+e) + l(ln(a2A2)+/), (3.6) 

•^LHQET = ('^■^) 

47iT = -^ + rf^'""^- (3.8) 

Using them, the matching factors that appear in Eq. (2.78) are 

-^CHQET - -^LHQET = ) + ^HQET - > + + , (3.9) 

"^CHQET "^LHQET ~ " ' {O.lU) 

l(ma) _ i{ma) _ _ i{ma) 
•^CHQET -^LHQET ~ " • \'^-^^) 

Infrared divergences cancel between CHQET and LHQET. We also note that the tree-level 
amplitudes of CHQET and LHQET are not generally the same. The finite values /, e, 
and S''"'^ are obtained fr om the following one-loop calculations. 

Light quark propagator correction 

The light quark self-energy terms from the rising-sun (RS) and tadpole (TP) diagrams are 




^r(f) = ' f * — ^ ' ^ ip) = ^ — . (3.12) 

P k + p P 




The wave function renormalization is 

^ g(sf(p) + sr(p)) 



2/=0 



\2 

ln-+/ 



(3.13) 



For the chiral light quarks employed in this paper there are no 0{a) contributions to this 
part. The numerical value of / for the domain-wall quarks is given in Refs. [31, 32]. We 
cite their results and quote them in Tab. 5. 
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Static heavy quark propagator correction 

The static self-energies from RS and TP contribution are represented by 




p k+p p 

and the radiative correction to the static mass is 




6M 



^hliPo = 0) - sTP(po = 0) ^ (±yCp6M. 



\AirJ 



, (3.14) 



(3.15) 



The wave function renormahzation is 
1 



d[j:^^^{po) + ^tlipo) 



1 + (1^) ' [-2 ln{a^X^) + e] . (3.16) 



Po=0 



As mentioned in Ref. [33], the static heavy quark propagator has no 0{a) error after 
imposing the on-sheU condition po = 0. The finite part e can be decomposed as 



e = n + 6M, 



(3.17) 



where TZ is defined in Eq. (B.15). 



Heavy- light vertex correction 

The heavy-light vertex correction has an 0(a) part which is determined by expanding in 
powers of the external momentum and the light quark masses. In this expansion, the equa- 
tions of motion are used. The expansion with respect to static quark external momentum 
always vanishes owing to the on-shell condition, whereas that with respect to light quark 
survives due to anisotropy between space and time in the system. The expansion of the 
vertex correction on the lattice has the form: 




k + p P 



(3.18) 



3aX 



+ 



in ■ p) 



+G 



4tt 
aX 



iTiq + 0{p'^,pmq,ml) \ + 0{g'^). 
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3.2 AB = 2 four-quark operator 

The one-loop transition amplitudes (2.81), (2.82) and (2.83), which are the relevant parts 
for the one-loop matching, are 

(f|OL|i)HQET = |l+ ((f|OL|i))HQET (3.19) 

+ (£)'i3(^'")((f|aO^^|i))HQET + (£)' ^("^'^)((f|aO;vAf |i»HQET + 0(/), 
(f|Os|i)HQET = ((f|05|i))HQET + 0{g^), (3.20) 

where B^^'^ = B^, B^P"^ = ^{^"^ and B^"'"'^ = respectively. Note that only the 

tree- level amplitude is needed for the amplitude of O5, since the coefficient Z2(/ib,/x) in 
Eq. (2.48) begins at O^g"^). In the expression above, the continuum values B^^\ QiP^) and 

4hqet = ^(cf + (" ^ + (^-^^^ 

iV, - 1 / „ A2 _ \ AT^ - 1 / _ a2 



+Cf (^-2 In ^ + EnQET^ + (^In ^ + F 

4'SU--(c. + #)(-|^). ' (3.22, 

4-,.-(c..^)(-i^), (3.23) 

where the values of L^hqet and F'hqet are shown in Eq. (2.47), F in Eq. (3.5) and also 

Chqet = 0, V = -5. (3.24) 

The lattice counterparts are 

€i1qet = ^{c, + ^) (- In (a^A^) + 

+ ^ (-'^^^ + ^) + ^ + 

+Cf (-2 In (a^A^) + e) + Cf (in (a^A^) + /) , (3.25) 

4SW-(^.+^)(-|^+^")- <^-^«> 

8ax^-(c..^)(-^^ + .'-'). (3.27) 
Using them, the matching factors that appear in Eq. (2.87) are 
4hqet - 4'i1qet = In (aV^) +2^Cf + (z^hqet - d^'^) (3.28) 

+ ^ (Chqet - c Fhqet - t") CiT' (-Ehqet - e + F - f) , 
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and (3.28), we need the following loop calcula- 
Light-light vertex correction 

For the AB = 2 four-quark operator, we need the light-light vertex correction: 

^r,, = =(T^>^^,^)(£)'[41n(a2A2)+^;]. (3.31) 

For domain-wall quarks, v in Eq. (3.31) was calculated in Ref. [32, 34]. We quote those 
results in Tab. 5. For domain- wall quarks there is no 0{a) error in this part even when the 
quarks are off-shell. 

Heavy- heavy vertex correction 

For the four-quark operator we also need the heavy-heavy vertex; 

h. = {T^jf:^T^7f:) [-2ln{a'X') + c] . (3.32) 

c = e-5M. (3.33) 
Note also that if we impose the on-shell condition, this part has no 0{a) error. 

4 Lattice perturbation theory 

The full one-loop lattice perturbation theory calculation is presented in Appendices A 
and B. Here we make some general comments on the calculation. 

4.1 Lattice action 

The lattice action comprises three pieces: 

•S" = S'static + 'S'dW + 'S'gluon, (4.1) 

where ^static is the static quark action describing the heavy (6) quark, 5dw is the domain- 
wall fermion action describing the light {u, d, s) quarks and S'giuon is the gluon action. This 
action satisfies the symmetries presented in Sec. 2.4.1. 



■"■CHQET "LHQET 



j2(ma) j3(-ma) 
"CHQET ~ "LHQET " 

To obtain c and v in the expressions (3.25) 
tions. 




k -k 

The finite part c is exactly 
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4.1.1 Standard static action with link smearing 

The standard static quark action [35] is given by 



5. 



static 



{h+{x) [h+ix) - U^ix - 6)h+{x - 6) 

X 

-h-{x) Uo{x)h^{x + 6) - h_{x) |. 



(4.2) 



The lattice derivatives used here are not symmetric in order to avoid fermion doublers. We 
introduce link smearing to improve the signal to noise ratio [5]. The modification is just 
to replace the link variable Uq{x) with a 3-step hypercubic blocked one Vq{x), defined by 



Vf,{x) = Projsu(3) 
V^,;uix) = Projsu(3) 
V^,;upix) = Projsu(3) 



1 - ai)[/^(x) + ^ V^,f,{x)Vf,,^{x + i))Vj.Jx + fi) 



(1 - a3)C/^(x) + ^ Y Urj{x)U^{x + ml{x + fl) 



±r]j^p,u,ii 



(4.3) 
,(4.4) 
(4.5) 



where Projgu(3) denotes a SU{3) projection and (01,02,03) are the hypercubic blocking 
parameters [36]. We use the parameter choices: 



(01,02,03) 



unsmeared (V^ = U^) 
APE with Q = 1 [37] 
HYPl [36] 
HYP2 [5]. 



(0.0,0.0,0.0) 
(1.0,0.0,0.0) 
(0.75,0.6,0.3) 
(1.0,1.0,0.5) 

4.1.2 Domain-wall fermion action 

The domain- wall fermion action with fermion mass parameter mj is ^ 

Ls _ 

s,s'=l x,y X 

DfJ^ix, y) = D\x, y)6ss' + D\s, s')5,y + (M5 - h)5,s'5^y. 
D\x,y) = Yl[il+lp)U,ix)6^+f.,y + il-j,)Ul{y)6., 



(4.6) 



x-ii,y 



PBh,s' {s = 1) 

Pb5s+i,s' + PlSs-i,s' {I < s < 
PlSls-i,s' (s = Ls) 



(4.7) 

(4.8) 
(4.9) 

(4.10) 



where ipsix) is a 4+1-dimensional fermion field. The fifth dimension extends from 1 to Ls 
and is labeled by the indices s and s' . The domain- wall height (fifth dimensional mass) 



^The expression of the action we employ here is different from that used in many simulations. See 
Appendix C. 
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Ms is a parameter of the theory which can be set between < < 2. The physical 
four-dimensional quark field q{x) is constructed from the 4-|-l-dimensional field ips{x) at 
s = 1 and Ls'. 

q{x) = PrM^) + Pl^lAx), (4.11) 
q{x) = Mx)Pl+^lSx)Pr. (4.12) 

In our perturbative calculations we assume that Ls is infinity such that the right and left- 
handed modes are decoupled and chiral symmetry is restored. We follow the formalism 
developed in Refs. [31, 38]. The physical quark propagator is written as 

S,{p) = {q{-p)q{p)) = "^'l (1 + 0{p',pmf,m})) , (4.13) 

where wq = 1 — M5. From Eq. (4.13) we see that the quark wave function has a domain-wall 
specific factor (1 — Wq)^/^ and that the relation between the quark mass ruq and the domain- 
wall fermion mass parameter mf is = (1 — WQ)mf. Hence the relationships between the 
lattice and continuum tree-level amplitudes in Eqs. (2.76), (2.78), (2.84), (2.85) and (2.87) 
are: 



5(J±r) = S{J±rD) = S{J±tm) = (1 - wl)-^'', (4.14) 
S{Ol) = S{Ond) = S{Onm) = S{Os) = S{Oj^) = 5(0^) = (1 - wl)-\ 

As is discussed in Refs. [31, 39], the domain- wall height M5 receives an additive quan- 
tum correction. This leads to a renormalization of the factor (1 — Wq), which we write as 
Zyj. The one-loop expression for Z^, is 

Zw = 'i- + 7T~v2^F~r~^^^ = 1 + -j^-t^Cfz^, (4.15) 
(47r)^ 1 - ""^0 (47r)^ 

where values for are listed in Ref. [32] and quoted in Tab. 4. This domain- wall specific 
factor can be treated as a part of the wave function renormalization. Note that the domain- 
wall fermion is automatically off-shell 0(a) improved [40]. 

4.1.3 Gluon action 

We consider a class of RG-improved gluon actions: 



5, 



gluon 



8ci) V ReTr[[/p] + ci V ReTr[[/R] ) , (4.16) 



where qq denotes the bare lattice coupling, Up is the path-ordered product of links around 
the 1x1 plaquette P and Ur is the path-ordered product of links around the 1x2 rectangle 
R. We will present numerical results for the following choices of ci: 

(standard plaquette action) 

-1/12 (Symanzik action) [41, 42] 
Cl = ■ (4-17) 

-0.331 (Iwasaki action) [43, 44] 

-1.40686 (DBW2 action) [45] 
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4.2 Mean field improvement 

Mean field (MF) improvement is important for precise estimation in lattice perturbation 
theory [46]. The lattice action is improved by replacing the gluon links with C/^/uq, 
where uq represents the MF value of U^. For the value of uq we choose the fourth root of 
the expectation value of the plaquette P taken from the simulation and to accomplish MF 
improvement we also need the one- loop perturbative expansion of uq: 

uo = P'/' = l-9'Cp^, (4.18) 
where the MF factor Tmf takes the values [32] 



TuF 



1/8 (standard plaquette action) 

0.0915657 (Symanzik action) 

0.0525664 (Iwasaki action) 

0.0191580 (DBW2 action) 



(4.19) 



First we present MF improvement for the light quark sector. For the domain-wall 
height M5 we have the following replacements [32]: 

Ms Mi''' = M5 - 4(1 - no), (4.20) 

woiM^) w^"" = woiMi"^) = wo{M5) + 4(1 - uq), (4.21) 

4,,,MF 

UM,) z^^{M,) = z^iMi"^) + \ [A^fTuF. (4.22) 

1 — (Wq Y 

Replacements are also required as in the usual (Wilson quark) treatment: 

/(M5) r\M,) = /(Mr) - (4vr)2^, (4.23) 

v{M,) 7;^^(M5) = viMi'''), (4.24) 

q uy\, (4.25) 



(1 - wl)mf mf^ = ^(1 - {wrr)mf, (4.26) 

Uq 



^. (4.27) 

MF improvement for static quarks was first discussed in Ref. [47] and is summarized 
here: 

(1) As mentioned in Ref. [35], the power divergent piece in the static heavy quark self- 
energy causes a change in the fitting function of the B meson correlator. Here, let 
us illustrate this issue using the static quark propagator. The propagator is given by 
Eq. (A.l) and its inverse is 

{Ship))-' = l-e-'P'>, (4.28) 
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where only the h+ part is taken. The renormahzation of the static propagator is 
calculated from the static self-energy E/j(po) (3.12) and the inverse bare propagator 
is written as 

(Ship))-' ~ Z~\l - e-'P^ + 5M) {po « 1), (4.29) 

where and 6M at one-loop are given by Eqs. (3.16) and (3.15), respectively. Then, 
the coordinate space static quark propagator is 

5/,(i) ~ ^he-i'^(^+^*^)(*+^) (t>l). (4.30) 

Here we consider two types of fit function for the propagator (4.30): 

Fi{t) = Ae-^^'+^\ (4.31) 
F2{t) = Ae~^K (4.32) 

Although the difference between them naively seems to be an irrelevant 0{a) effect 
in the continuum limit, it requires special attention owing to the 1/a divergence in 
6M. When we choose the fit function (4.32), Eq. (4.30) is rewritten as 

where = Z}i/{1 + (5M), which leads to a "reduced value" for e (3.17): 

eR = e- 5M. (4.34) 

We use the fit function (4.32) with the reduced value en rather than (4.31) to avoid an 
unwanted 0{a) error in the light quark propagator which obeys the function (4.32). 

(2) The link variables in the static quark action are smeared. One way to implement MF 
improvement for them is to replace the smeared gauge links Vq with Vq/vq, where vq 
represents an MF value of Vq [48]. Associated with this, the one-loop perturbative 
expansion of vq is introduced using the expectation value of the plaquette with the 
link smearing Psmeared: 

T-ismeared 

= (^'smeared)'/' = 1 " Q^Cp^^ . (4.35) 

(3) In coordinate space, the static quark propagator Sh{t) is related to the MF-improved 
propagator Sh{t) by 

Sh{t) = vlSh{t), (4.36) 

since each link variable in the propagator is multiplied by 1/vq in MF-improved 
case. The relation (4.36) is satisfied nonperturbatively but we show the perturbative 
derivation for pedagogical reasons. Starting from the static quark action without MF 
improvement, the action with MF improvement is obtained by rescaling the static 
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1 /2 

quark field by Vq and adding a static quark mass Mq = {1/vq — 1). The interesting 
feature of the static quark is that the mass can be absorbed into a shift of the external 
momentum p in the propagator, whose momentum is k+p {k is the loop momentum), 

{S^h (k + p; Mo)r' = 1 - e^^(^«+P«) + Mo 

= {l + Mo){S±h{k+p';Mo = 0)y' (4.37) 

where Po = PoT^ln(l + Mo), a shifted external momentum. Perturbative calculations 
are not affected by this mass shift but the position of the pole in the propagator is 
changed. The value of e is not altered but the static quark mass correction changes: 

6M{Mo) = (1 + Mo)6M{Mo = 0). (4.38) 

Thus we find the same result as Eq. (4.36): 

S,it) ^ ^Z,e-'<'^^'^<^o-')>'' = y^,S,{t). (4.39) 

Combining the results above, we relate the coordinate space quark propagator to the MF- 
improved one by 

^,^_ln(l+5M)i ^ ^,g-ln^(l+5M)i^ 

in which Zf^ and 5M are estimated using the MF improved coupling guF- The replacements 
needed for MF improvement are 

c^c^^ = c, (4.41) 

eR eT = en, (4.43) 
h — > h, (4.44) 

1 f , 2 / T^smeared \ ^ 

m SM"^^ = JL J 1 _ „o + -^Cf Um - (47r)2iML j I .(4.45) 

At one-loop in MF-improved perturbation theory, the renormalized coupling in the 
continuum MS scheme g^ifJ-) at scale fi and the bare lattice coupling go are related by 

IP 22 / 4 \ 

^ = ^ + <i, + c, + ^ H,a) + [df - ^ H,a) j , (4.46) 

where dg and Cp are dependent only on the gluon action and dj is dependent only on the 
fermion action. Cp is obtained from 

Cp = 2CfTmf, (4.47) 
and the value of dg is given by [42, 49-54]: 

-0.4682 (standard plaquette action) 

-0.2361 (Symanzik action) 

0.1053 (Iwasaki action) 

0.5317 (DBW2 action). 



dg 
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Values df for domain- wall fermions are given in Ref. [55] and quoted in Tab. 6. 

Collecting the pieces needed for the MF improvement, we now have final expressions: 



j: 



CHQET 



±r 



(l-K^^)2)-V2(^MF)-l/2^1/2 



(4.49) 



y{l)MF ^LHQET ^ 



(pa)MF 



tLhqet , fr xLhqet 



Uq 



tLHQET (pa) .LHQET ^(^a) ^HQET 
^_I_P -r u<v_|_p£) -r cp ",<^_|_p^ 



CHQET 



(1 - K^^)2)-n^r)-^-o 



(4.50) 



^(l)MF^LHQET _^ 



,(pa)MF 



no 



^LHQET 



,(ma)MF 



^LHQET (pa) ^LHQET ^(ma)^ 



Uq 
LHQET 



LHQET 
NM 



o 



CHQET 



(l-«-)^)-(zr)-^noO^ 
where J^p^^'^ and O^^^'^ are introduced by 



(4.51) 



J; 



LHQET 



itEAf 



LHQET 

NM 



MF^2 



7LHQET 



J. 



o 



LHQET 



In the above equations the matching coefficients are written as 



7(1)MF 

^r(i?) 

^(pa)MF 
Zp 

7 (ma) MF 



\47r 



Cf 



-^CHQET — -^LHQET 



(1)MF 



JL)- 

A-k) 
An J 



Cf 



Cf 



-^CHQET — -^LHQET 
-^CHQET — ^LHQET 



(pa) MF 
(ma)MF 



\47r 



(4.52) 



CfzP^^"", (4.53) 



(^Ycfz. 



(pa)MF 



V47r/ 
V47r/ 



Cfz- 



(ma)MF 



z 



(1)MF 
L{R) 



V47r 



Z 



7(pa)MF 



(ma)MF 



47r/ 
An J 



^CHQET — ^LHQET 

-I R 

2 r . . ] (pa)MF 

^CHQET — ^LHQET 



(1)MF 



1 ^ f^V 9^^^ 



V47r 



.(pa)MF 



^CHQET — ^LHQET 



(ma)MF 



g Y .(ma)MF 

4^) 



(4.54) 
(4.55) 
(4.56) 
(4.57) 
(4.58) 



where g is replaced by (/^5g(/^ = a ^) in Eq. (4.46), and is used instead of e in the 
matching factors (4.53) and (4.56). The numerical values of the one-loop coefficients in 
Eqs. (4.53)"(4.58) are shown in Tabs. 24 - 35. 

4.3 Alternative form for the 0{pa) terms 

We defined 0{pa) operators containing covariant derivatives for the quark bilinears and for 
the AB = 2 four-quark operator in Eq. (2.21) and Eq. (2.28) respectively. These operators 



-24- 



can be conveniently rewritten at 0{g^) using the equations of motion as: 



J±rD = (G(=Fno5o + I - uq) - mg)J±r, 



Ond = 2uo[h+^l:q] [do - do) [h.j^q] - 2(1 - uo)Ol - Onm, 
where we used the tree-level MF-improved equations of motion for the light quark 



Uo Uo 



and for the static quark 



h± 



nsmcarcd 
n 



0, 



1^1 

vo 



0, 



(4.59) 
(4.60) 



(4.61) 



(4.62) 



(The unimproved case is obtained by setting = 1 and = 1.) In the derivation of 
Eqs. (4.59) and (4.60), we replace vo with uo because the difference between the smeared 
and unsmeared gauge field is O(a^) [56, 57] and it is irrelevant for the 0{a) improvement 
program. We also note that there are derivatives rather than covariant derivatives in 
Eqs. (4.59) and (4.60). The cancellation of gluon fields occurs between the covariant 
derivative of the static and that of the light quark: 



smeared , 
J^O + 



i^^— d ^ + [do + igAo + 0{g'a) 

ig{l + 0{a^))Ao + 0{g^a)) + (^o + igAo + 0{g^ 



do + do + 0{g a,ga ). 



(4.63) 



This cancellation is valid at full order in g and the 0{a) deficit in Eq. (4.63) becomes an 
O(a^) error in the 0(a) improvement program. 

Now we apply the operator J±rD (4.59) to the B meson correlator: 



(J_rz)(t)(J-r(0))n 



G 



-uoE^'""^ + 1 



no) -m,}(J_r(t)(J-r(0))^), (4.64) 



where E^^'^'^ denotes the B meson binding energy, which comes from fitting the correlator 
to a function proportional to e"^*"'"**. The total time derivative in J-rn picks up -Ebind 
which contains a 1/a divergence originating from the static quark mass shift 5M. At first 
glance, this appears to conflict with using J±rD, the original 0{pa) operator, because the 
spatial derivative in J±rD seems not to pick up the 1/a divergence in the static quark 
momentum. We will discuss this point in more detail in next subsection. While Eq. (4.64) 
is a part of 0{a) term, a combination of the 0{a) and the 1/a affects the 0(1) term. This 
is, however, an effect at higher order in the one-loop perturbative calculation. 

4.4 Power divergence from 0{pa) operator 

As already mentioned, we need to pay special attention to the 1/a power divergence when 
taking the continuum limit. One issue is the fit function ambiguity of the correlator, 
which we discussed in Sec. 4.2. Another issue is the power divergence coming from the 
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0{pa) operator, mentioned above in Sec. 4.3. The loss of Lorentz invariance in the lattice 
discretization allows the 0{pa) operators to give rise to a power divergence. The divergence 
comes from quantum corrections, that is, 0{g^) effects, but because these operators appear 
at 0{g'^) in the 0(a) improvement terms, the power divergent piece affects only higher 
orders in perturbation theory. It is, however, useful to see how this divergence arises 
from the 0{pa) operators both for the higher order perturbative and for non-perturbative 
renormalization. 

Consider first the 0{pa) operator for the quark bilinear J±rD in Eq- (2.21). To locate 
the power divergences, we calculate the transition amplitude of this operator at 0{g'^). 
When the on-shell condition is imposed, the external quark momenta should be the renor- 
malized ones. The external light quark momentum does not cause the power divergence, 
since the massless point is well defined for chiral fermions. Note that even for the Wilson 
fermion action, which has intrinsic chiral symmetry breaking, a power divergence can be 
subtracted by adjusting the massless point using the critical hopping parameter Kc- The 
on-shell external static quark momentum p' , however, receives a static quark mass renor- 
malization, leading to ±«Po + SiM = 0, which means p' includes a power divergence. In the 
analysis up to one-loop we include tree and one-loop diagrams. Some tree diagrams acquire 
a power divergence in the on-shell external momentum p' and these tree diagrams need to 
be taken into account. Divergences from p' are also contained in the one-loop diagrams, 
but are a higher order effect by and therefore can be omitted in our analysis at the 
one-loop level. 

For the operator J±td the tree- level amplitude does not pick up a power divergence 
since the operator contains a derivative acting only on the light quark: 




i [0 + 0{pa)] . 



(4.65) 



h±T(f ■ D)q " 

The following one-loop diagrams need to be calculated: 




(4.66) 




(4.67) 




(4.68) 
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The power divergent piece comes from the loop diagram of Eq. (4.66). Again, this power 
divergence contributes at 0{g^) which is higher order than the one-loop calculation in this 
paper, since the 0{pa) term is already 0{g'^). 

In Sec. 4.3, the transformation of the original 0{pa) operator to a total time deriva- 
tive operator was discussed. Using the equations of motion the operator ^do{h±T^oq) is 
obtained. This operator gives a power divergence from the tree diagrams with the on-shell 
condition: 



Tdo{h±T'yoq) 




-[-G5M + 0{g\pa)], (4.70) 
a 



= -[0 + Oig\pa)]. (4.71) 

The loop diagram (4.71), however, does not produce a divergence owing to the nature of 
total derivative at 0{g^). The power divergence is the same as that in the original operator 
J±rD which is consistent with the fact that the equations of motion do not contain any 
power divergence. (The power divergent structure in the equations of motion is presented 
in Appendix D.) Now we answer the puzzle introduced in the last paragraph in Sec. 4.3. 
The power divergence in the operator T<9o(^±r7o?) comes from the on-shell static quark 
momentum, while the divergence for the original 0{pa) operator J±rD arises through 
quantum correction to its vertex function. 



5 Numerical value of the matching factor: an example 

In this section we show an example how to obtain the numerical values of the renormal- 
ization constants and 0{a) improvement coefficients using the simulation settings from 
Ref. [58]. The simulation parameters are presented in Tab. 1. 

5.1 Continuum matching factor 

We fix the coupling constant at the Z-boson mass to the PDG value as{mz = 91.1876 GeV) = 
0.1184 [15], from which Us at different scales are calculated using four-loop RG run- 
ning [59, 60]. Because the simulation is performed not including dynamical charm quark, we 
employ a two step running as mentioned in Sec. 2.5. The obtained coupling constants are: 
a,(mfe = 4.19 GeV[15]) = 0.2260, a^(mc = 1.27 GeV[15]) = 0.3919, 0,(0"^ = 1.73 GeV) = 
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Lonngurations (Aj 


Coniigurations [d) 




2.13 


2.25 


T,^ X T X T, 


X fid- X Ifi 


39'^ fi4 X Ifi 


(GeV) 


1.73(3) 


2.28(3) 




0.003152(43) 


0.0006664(76) 


Ms 


1.8 


1.8 


P (chiral limit) 


0.5883 


0.6156 




1.3032 


1.3432 



Table 1. Simulation data for 2 + 1 flavor dynamical domain- wall fermions with Iwasaki gluons 
from the RBC-UKQCD CoUaborations [58]. 



0.3204 and as{a ^ = 2.28 GeV) = 0.2773. Using these couplings, the matching factors 
between CQCD and CHQET, (2.53) and (2.54), are: 







Cr=7o7s 


(mb,a 


-1 = 1.73) = 






Cr=7o75 


(mfe,a 


-1 = 2.28) = 








T 




T r 


Zi{mb 


= 1.73) 




0.7483 




Z2{mb 


a-i 


= 1.73) 




-0.1439 










T 




T r 




a-i 


= 2.28) 




0.7483 




Z2{mb 




= 2.28) 




-0.1439 





0.9520 X 1.1550 x 0.9521 = 1.0470, 
0.9520 X 1.1550 x 0.9196 = 1.0112, 



(5.1) 
(5.2) 



1.3345 
-0.0526 1.0921 

1.3345 
-0.0526 1.0921 



0.9055 
0.0141 0.9706 

0.8442 
0.0231 0.9500 





0.9088 




-0.1525 




( 




0.8457 




-0.1493 



(5.4) 

5.2 Lattice to continuum matching factor and 0{a) improvement coefficient 
in HQET 

We present one-loop perturbative results for the lattice to continuum matching factor and 
0(a) improvement coefficients in HQET. The expectation value of plaquette P is the value 
in the chiral limit, obtained by linear extrapolation in quark masses using the lightest two 
degenerate up and down quark mass parameters. For Configurations (B) the deviation 
from that using three quark mass parameters is less than 0.01%. 

Given the value of P, the matching factor and 0{a) improvement coefficients are 
obtained by the following steps: 

(1) The MF value of link uq = P^/^ and the MF improved value of the domain- wall height 
Mg^^, which is calculated from Eq. (4.20), are obtained as presented in Tab. 1. 



(2) The value of Ti^ at M5 

MF 



Mf^ is taken fr om Tab. 4 using cubic spline interpolation. 



Then z^''''^ (M5) is obtained using Eq. (4.22). The values are presented in Tab. 2. 
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5ms /47r 

{1)MF 
T(R) 
(pa)MF 



-(ma)MF 

41)MF 

^L(R) 

.(pa)MF 

~(ma)MF 



2:, 



MF 



{1)MF 
V(R) 
-v{pa)MF 



2;, 



^1 



(7na)MF 



Z 



{1)MF 
L(R) 
.(pa)MF 

(ma)MF 



unsmear 



Configurations (A) 
APE HYPl 



HYP2 



-4.87061(56) 
0.41612(41) 
1.70169(57) 
-15.4170(15) 
-0.83223(82) 
-3.4034(11) 



0.1769 
5.257 

-1.58382(56) -1.50401(56) 



-3.48036(41) 
-4.12249(57) 
-4.4611(15) 
6.96072(83) 
8.2450(11) 



-3.68664(42) 

-4.43848(57) 

-4.1950(15) 

7.37328(84) 

8.8770(11) 



0.07768(56) 
-6.41186(44) 
-8.52774(57) 

1.0773(15) 
12.82372(87) 

17.0555(11) 



unsmear 



Configurations (B) 
APE HYPl 



HYP2 



-4.88040(68) 
0.41254(50) 
2.13243(68) 
-15.3875(18) 
-0.8251(10) 
-4.2649(13) 



0.1683 
6.133 

-1.58901(68) -1.50890(68) 



-3.48553(51) 
-3.84087(68) 
-4.4162(18) 
6.9711(10) 
7.6817(14) 



-3.69280(51) 
-4.15412(68) 
-4.1492(18) 
7.3856(10) 
8.3082(14) 



0.07608(68) 
-6.42056(53) 
-8.33056(68) 
1.1340(18) 
12.8411(11) 
16.6611(14) 



Table 2. Numerical values of the one-loop coefficients of matching factors (7£jg/47r, z. 

(ma)MF 



-(pa)MF ^(ma)MF «(1)MF ^{pa)MF 



jMF 5(1)MF 



and 



(3) The one-loop coefficient of the matching factor (without the renormalization factor of 
M5) and 0(a) improvement coefficients are read from Tabs. 30-32 using cubic spline 
interpolation. The values are presented in Tab. 2. 

(4) The coupling (7j^(/i=a~^) is obtained from Eq. (4.46) as presented in Tab. 2. In the 
calculation the value of at M5 = is used, which is read from Tab. 6 using 
cubic spline interpolation. 

(5) The matching factors and 0(a) improvement coefficients in Eqs. (4.49)-(4.51) are 
obtained by combining the steps above and are presented in Tab. 3. 

As expected, for the 0{aP) matching factors Z^^^ and -^^j^^ link smearing tends to 
produce a small correction; HYP2 smearing gives an especially tiny change. For the values 
of the one- loop 0{pa) and 0{ma) improvement coefficients Z^'^^^^ , ^{pa)MF ^ ^^ma)MF 
2,{rna)M¥ ^ link Smearing, especially HYP2 smearing, tends to produce large coefficients. 
This property depends on the choice of the gluon action, however, since when we choose 
the standard plaquette action, the link smearings tend to give tiny coefficients. 
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Configurations (A) 




Configurations (B) 






unsmear 


APE 


HYPl 


HYP2 


unsmear 


APE 


HYPl 


HYP2 


-L jU JO 


0.8513 


0.9091 


0.9105 


0.9383 


0.8684 


0.9243 


0.9256 


0.9526 


(pa) 

, '-P " 

{ma) 

Zl 


-0.0089 
-0.0331 
0.6873 


0.0746 
0.0802 
0.8227 


0.0790 
0.0864 
0.8260 


0.1374 
0.1660 
0.8911 


-0.0083 
-0.0379 
0.7184 


0.0703 
0.0683 
0.8514 


0.0744 
0.0739 
0.8546 


0.1294 
0.1482 
0.9187 




-0.0134 


0.1119 


0.1185 


0.2061 


-0.0125 


0.1054 


0.1117 


0.1942 


(ma) 


-0.0497 


0.1203 


0.1296 


0.2489 


-0.0569 


0.1025 


0.1108 


0.2222 


Zs 




0.9645 






1.0040 





Table 3. Numerical values of the one- loop lattice to continuum matching factors and 0{a) 
improvement coefficients in the HQET (MF-improved value). 



5.3 0{a) corrections for physical quantities 

In this subsection we estimate the 0{a) corrections for B meson quantities using the results 
obtained in Sec. 5.2. Using (4.59) the \B) to |0) transition amplitude leads 



/nl /iO(a)impr| 

(01^0 |i^)cHQET ^ >a) / ^bind ^ ^ W / M _ >a) \ 

(0|Ao|S)cHQET r=7075 V " B t U y ^ r=7o75 r=7o75; / 

(5.5) 



where E^^'^ is the binding energy of the B meson obtained from the fit of the correlator 
(^o(i)-4j(0)). Note that one does not need to compute (A[f"V)Aj(0)). For the AS = 2 
four-quark operator, we need to calculate the matrix elements of the operators oj^"^ and 
0^"\ For simplicity we use the vacuum saturation approximation (VSA), where the \B) 
to \B) transition amplitude can be written in the form: 

/^m'^(")™p''i R\vsA 

I-"/CHQET _ (pa) / pbind , „ -, \ , 9 ( (ma) (pa) \ 

(5.6) 

The fact that the overall matching factor in Eqs. (5.6) is the square of that in Eq. (5.5) 
is consistent with the spirit of the VSA. Using Eqs. (5.5) and (5.6), The 0{a) corrections 
to the B meson decay constant and /\B = 2 matrix element are evaluated. The 0{a) 
correction to the B meson decay constant is given by Eq. (5.5). For the Ai? = 2 matrix 
element Mb-, we can estimate the 0(a) correction using Eq. (5.6), where, of course, the 
bag parameter Bb oc MbIb'^ VSA has no 0{a) correction. The 0{a) correction to 

the SU{3) breaking ratio ^ = fBs\/ BbsI {fB^^J Bb^) in the VSA is 

'^^) = 1 + 4it.a-o(i^r - ) + - cSU,) (-,(.) - m,(d)). 

? /VSA ^ ^ 

(5.7) 
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We apply the above equations to the actual numerical simulations performed in Ref. [29]. 
The (3 in the simulation is the same as for configurations (A) in Tab. 1 but the lat- 
tice size is X r X Ls = 16'^ x 32 x 16. The B meson binding energy is roughly 
£;bind ^ Q_g (APE), 0.5 (HYP2). To estimate the 0{a) contribution, we omit the 0{ma) 
effect because the domain- wall mass parameter nif in the simulation is 0(0.01), which 
is much smaller than that of The conclusion is that the 0{g'^a) effects on fs are 

around 3%(APE), 4%(HYP2), and on Mb are around 6%(APE), 8%(HYP2) in the VSA. 
Using the assumption — -E'^™'^) ~ {itiBs ~ ''^B^)^ the effect on ^ is less than 2% in 

the VSA. 

6 Conclusion 

In this paper, we have calculated renormalization factors for the heavy-light quark bilinear 
operator and AB = 2 four quark operator in the static heavy and domain-wall light quark 
system, including 0(a) corrections. Even for domain-wall fermions, which have good chiral 
symmetry, there is an 0{a) correction. The allowed set of 0{pa) and 0{ma) operators was 
constrained by symmetries. We showed the one-loop perturbative calculation of the lattice 
to continuum matching factors and the 0(a) improvement coefficients taking into account 
link smearing (APE, HYPl, HYP2) in the static heavy action and four types of gluon action 
(Plaquette, Symanzik, Iwasaki, DBW2). The results showed that the 0{pa) correction is 
not negligible and should be included in simulations used for precise determinations of 
CKM matrix elements. 
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A Lattice Feynman rules 

In this appendix we list the Feynman rules used in the lattice perturbation theory. 

Static heavy quark sector 

• static/anti-static heavy quark (+/— ) propagator 

Sh±{k) = {h±{-k)h±{k)) = ^—i^Pi, (A.l) 
where P± = (1 ±7o)/2. 
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• static/anti-static heavy quark (+/— ) ~ one gluon vertex 

W^^{k,k') = ^igT^5o^e^'^^°+^'^^/'^P±. (A.2) 

• static /anti-static heavy quark (+/— ) - two gluon vertex 

• hnk smearing 

For the static heavy action, we impose Hnk smearing. This smearing changes the 
Feynman rules [57, 61]. In this case the rules of the static heavy quark - gluon vertex 
(A.2) and (A. 3) are changed as 

W^^{Kk') W^\k,k') = Y,Kp{mp''{hk') = h^o{k)W^^{k,k'), (A.4) 

p 

W^^^'iK k') W^J'^'ik, k') = Kp{k)W%^^{k, k')Ku{k') (A.5) 

per 

= Ko{k)W^^'^''{k,k')hov{k')^ 



where 



h^uik) — 6^1, 



+ ^kf^Kn^^^ik), (A.6) 



p 

n^,{k) = l + a2{l + as)-^{l + 2a3){P-kl-kl) + ^ J] ^^'^^ 



v¥^p.,v 



kfi 



2sin(^), P = Ykl, (A.^ 



with smearing parameters {ai, 02,013). 

Domain-vifall light quark sector 
• light quark propagator 

Sg{p;mf) = {q{-p)q{p)) (A.9) 
—27^ sin + {I — VF(p)e~°(*'))mj 



z |<+ (1 - w'^)mf 
where 



(1 - Ty(p)e"{p)) + mjil - VF(p)e-"(p)) 
1 + 0{p'^,pmf,mj)) 



wo = l-M5, (A.IO) 

W{p) = wo + ^{1- cos p^), (A.ll) 
p 

W{p) cosh a{p) = i ^1 + VF2(p) + ^ sin^ p^ . (A. 12) 
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light quark to domain-wall fermion propagator 

{q{-p)i^s{p)) = 5,(p;my)(e-"(P)(^=-^)PR + e-°(^')(^-i)PL) (A.13) 

(M-pMp)) = (e-"(^=^^)PL + e-"(^-i)PR)5,(p; mj) (A.14) 
^(g-a(L.-.)p^ ^ e-«(^-i)p^)e--(p) (^Sgip;mf)mf - 1) . 

external light quark to domain-wall fermion propagator 



(gext(-p)V's(p)) 



+ {1 — WQ)mf 



{i^s(.-p)qc^t{p)) = Hs{p;mf) 



l-wi 



+ (1 — WQ)mf ' 



(A.15) 
(A.16) 



where 



Hs{p;mf) (A.17) 

[w^'-'Pr + t/^o~^^i) + rufWoiw^'-'PL + -u^r^^^?) + 0{p'^,pmf, m}) 
Hs{p;mf) 



(A.18) 



in which equation of motion for light quark is imposed, 
domain-wall fermions - one gluon vertex 



y^ij>,p')s 



■ , {p + p')^ , • ■ 

-igl I 7^ cos — ^+«sm 



A I yP^ P , ■ ■ {P + P')i^ \ . 



domain-wall fermions - two gluons vertex 



1 



sm ■ 



{P + P')^ 



+ cos 



{p+p'), 



some useful definitions 

S^{p;mj) = Y^{q{-p)lp^(p)){wL.-^PR + W^-'PL) 

s=l 

= -5g(p;m/)5^(p;0)e"(P), 

Ls _ 

S^{p;mf) = Y,{Qi-p)i^sip))iw^'~'PL + w'Q-'PR) 

s=l 

= {1 - Sq{p;mf)mf)Si,{p;0), 



(A.19) 



(A.20) 



(A.21) 



(A.22) 



where 



Sw{p; 0) 



Wo — e 



a{p) • 



(A.23) 



-33- 



Improved gluon sector 

In our calculation, the RG-improved gluon action is used. For the perturbation theory, we 
fix the gauge to Feynman gauge = 1). 

• RG-improved gluon propagator [42, 51] 

1 



where 



l-S. 



(A.24) 



Aik) 



(A.25) 



(A.26) 



p - c, U^f + E + 1 ((^')' + 2 E - E 

\ T J \ T T / 

-4c?E^pn^'- 

P ^7^P 

• some useful definitions for the smeared link 

p p 

^fuA;)^E^pp(^)^p-(^)^-(^)- 

p,(T 

B Calculation of lattice perturbation 

Static heavy quark propagator correction (Eq. (3.14)) 

The RS and TP contributions to the static heavy quark self-energy can be written in the 
one-loop integrals: 



(A.27) 
(A.28) 



S^^(p) = / D';,l{k)W^''{p,p + k)S±h{p + k)W^''{p + k,p), (B.l) 

J k 



T^llip) = \ jD%{k)Wt^^\p,p), 



(B.2) 
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where = J^^ (f^' order to obtain the radiative correction of the static mass and 
wave function renormahzation, we expand them hke 



, 2 (47r)2 



(B.3) 



Po=0 



^iiip) = 5]^r(po = 0) + 



-p± 



TP/ 



_9_ 

47r 



2 V ^ ^ 

2 



2 (4^)2 - 



(gTw _ 1) + 0((e^*P° - 1)') 



(B.4) 



Po=0 



-P± C-F^Tf^ - P± (£) C^il^TfHe^*^° - 1) + ©((e^^f" - l)^), 



where we define integrations: 

jin ^ -jz^i-P^ + ^PL-div ^ j^dW^ 



Jk L 



-ifco 



7^ 



PL-div 



(4vr) 



1 



(1 - e-^'^o +ae)2' 



(fc2 + A2)3/2(4 + fc2 + ^2)1/2 2(fc2 + A2)3/2 



R 



div 



(47r)2 /■ 0(1 -fc2) 



2 A (fc2 + A2)3/2 



4 (In 2 -In A - 1) 



rpCl 

-'4 



ncl-PL 



A: 

ici-PL I rpPL 

3 + ^3 ' 



D^},iO,k)-D^l-^iO,k) 



T. 



PL 



1 



where = • The numerical values of TiF^ and T^^ are 



7^ 



PL-div 



5.7531708(67), T^^ = 0.2527296(13) 



(B.5) 
(B.6) 

(B.7) 

(B.8) 

(B.9) 

(B.IO) 
(B.ll) 

(B.12) 
(B.13) 



in which their numerical numbers was also given in Ref. [12]. The new ingredients regard 
with regard to these quantities are the calculations taking into account the link smearing. 
The numerical values of and T^^~^^ are listed in Tabs. 7 and 8. Using them, the 

radiative correction to the static quark mass (3.15) and the wave function renormalization 
(3.16) are 



5M 



(47r)2 



\47r/ 



and then TZ in Eq. (3.17), which is a finite part of R"^^ , is given by 

n = R^^ +AlnX. 



(B.14) 
(B.15) 
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Heavy-light vertex correction (Eq. (3.18)) 

The one-loop vertex correction including 0{pa) and 0{ma) parts is calculated as 

5Th±q{p]mq) (B.16) 



s=l 



jk + 2p), 
2 

c + 2p)^ 



X {-i)gT I (5-;^(p + A;; m/) + + A;; 0)'u;om/) 7^ cos 
+ {Sw{p + k;mf) + Sy.{p + k; 0)womf) i sin - 



+ 0{p^ ,pmf,m?f) 



,pmf,mf), 

where we decomposed the vertex correction into 0{1), 0{pa) and 0{ma) parts in the last 
line. These 0{pa) and 0{ma) parts are the new calculations in this work, as well as the 
inclusion of link smearings. The 0(1) part is 



6r 



(1) 

hq 



•y-ci— PL T-PL— div rdiv •v-PL 



with integrations 



7-5 -PL 



(47r) 



kWo- e^C^) W{k) - e-^e^) ko 



(B.17) 
(B.18) 



rPL-div 



j;[z)S>)-<L(fc) 
1 



/cp cos 



2 / kp 



div 
X 

PL 



(4vr)2 
-(4vr)2 



^(l-fc2) 



1 2^feo\ 



fc k'^{k^ + {aXy) 
i^^f\j^SMD^o=\k). 



(B.19) 
(B.20) 
(B.21) 



The numerical values of and are listed in Tab. 11, and also that of X'^ is 



listed in Tabs. 12, 15, 18 and 21. The 0{pa) part is 
with integrals 



hq 



/jci— PL /yPL— div Tdiv ^ci— PL /jPL 



^ci-PL 
^X 



(4^)2 



6 



2 cos 



2 / 



+ 



5^(fc;0)e"W 
_ 1 - T^(fc)e"(*') 



(B.22) 
(B.23) 



1 - W^(fc)e"W 



sin^ /cp 



Pt^O 



+ 



^^-^ + 5^(fc; 0)e"('^A y sin^ kp 

1 - M^(/2)e°W ""^ ' > y ^ ^' 



2 VF(fc) + cos /cp — cosh Q(fc) 



W{k) sinha(fc) 
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J. 



PL-div 
X 



(4vr) 



+ 



1 S^(fc;0)e"W 
- ty(fc)e"{'=) 

1 



1 - VF(fc)e°(*') 
4^(1 - fc2 



^ gin^ A; ^^^^ ~^ ^Z' ~ ^^^^ '^^^^ 



E 



VF(fc) smha(fc) 



(fc2)2 



7' 



rdiv 
cl-PL 



PL 



-(4vr)2 

{±? 
2 

(47r)2 



Svr 



16 



Syj{k', 0), 



(B.24) 

(B.25) 
(B.26) 
(B.27) 



The numerical values of J'^^ and JT^'^ are listed in Tab. 11 and the those of J'^^ 
and J^"!"^^ are listed in Tabs. 12, 13 15, 16 18, 19, 21 and 22. The 0{ma) part is 

^(ma) _ _ ^ci-PL _ ^PL-div _ rdiv _ ^ci-Ph _ ^PL 



5T 



hq 



rrCx- 



(B.28) 



"PL— div JV'l'^ 

X ' 



ci-PL 



with integrals 

(4vr)2 
2 

2 

2 

(4vr)2 



^X' 



^PL-div 
^X' 



T-ai 



div 



' 1 5^(fc;0)e°W 
fc2 1 _ VF(fc)e"W 

0(1 _ 2^ 



1 - VF(fc)e"W ^ 
0(1 - fc2) 



(fc 



4, 



2^2 



-PL 



1 + 2S'^(fc;0)e 



^(^1^-^ 1 — VK(fc) cosh a{k) 



1 - VF(fc)e"W 



(B.29) 

(B.30) 

(B.31) 
1 



/c 



PL-div 



(4vr) 



P<|l + 25^(fc;0)e°(*') 



1 - VF(fc)coshQ(fc) 



(B.32) 
0(1 -fc2) 



1 - VF(fc)e"('^) J 1 - u^o (^^) 



div 



(4vr) 



0(1 - fc2 



2 Vfc fc2(/,2 + 



aX 



(B.33) 
(B.34) 



The numerical values of J^,^~'^" and /CP^-^^'^ are listed in Tab. 11 and the those of J^}~^^ 
and /C^i-PL are listed in Tabs. 13, 14 16, 17 19, 20, 22 and 23. Note that the first five terms 
in the right hand side of Eq. (B.28) come from the 0{pa) part by applying the light quark 
equation of motion —i^opo = i-y ■ p + niq in the calculation of the transition amplitudes. 
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C Relations of the Feynman rules and the one-loop amplitudes between 
different domain-wall fermion actions 



In this paper, we calculate the one-loop matching factors using the domain-wall fermions 
action defined in Sec. 4.1.2. In the most of numerical simulations, however, the following 
different action is employed : 



5'dw = ^ '^'4}s{x)D^^ {x,y)il)s'{y) -'^mfq{x)q{x), 

s,s'=l x,y X 

D^J{x,y) = D\x,y)5,s'+D\s,s')5^y + (M5 - b)5,s'6^y. 
D\x,y) = ^1 [(l-^^)C/^(x)5,+^,, + (l + 7^)f/t(y)5. 



-'x~ti,y 



[ PlS2,s' {s = 1) 

q{x) = PlMx) + Pr^^lM, 
q{x) = i^^ix)PR+lPrix)PL- 



(C.l) 

(C.2) 
(C.3) 

(C.4) 

(C.5) 
(C.6) 



When nif = +1, the anti-periodic boundary condition in the fifth direction, s-direction, 
is imposed for this action while the action in Sec. 4.1.2 imposes the periodic boundary 
condition. The Feynman rules for this action are different from those in Appendix A (The 
new rules from the action in this appendix are specified by prime symbols.): 



light quark propagator 

Sq{p;mf) = {q{-p)q{p)y = -Sq{p;mf). 

light quark to domain-wall fermion propagator 



(C.7) 



{q{-p)^sip)y = -5,(p;m^)(e-"(^')(^^-^)PL + e-°(^')(^-i)Pfi) (C.8) 
+ {Sg{p;mf)mf - 1) e-"(?')(e-"(^')(^^-^)PR + 6-""^^^^'-^^ Pl), 

(M-pMp))' = -(e-"(^^-^)Pfi + e-"(^-i)Pi)5,(p;mj) (C.9) 
^^^-a(Ls-s)p^ + e-°("-i)PR)e-"(?') {Sqip;mf)mf - 1) . 



external light quark to domain-wall fermion propagator 



(goxt(-p)V'.(p)>' 
(V'.(-p)^cxt(p))' 



1 - U^n ^/ 



-|- (1 — WQ)mf 



1 — w, 







if^ + (1 — wr\)m 



(C.IO) 
(C.ll) 
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where 



(C.12) 



m. 



H's{p;mf) (C.13) 
{w^'-'Pr + w'q-^Pl) - {w^'~'Pl + w'^-^PR)womf + 0{p'^,pmf, m)] 



domain-wall fermions - one gluon vertex 



<(p,A,s' = -i5T^(-7M 



, cos ■ 



domain-wall fermions - two gluons vertex 



{p + p'] 



+ 1 sm ■ 



5s,s'. (C.14) 



sm ■ 



{P + P')^l , {P + P') 



+ cos ■ 



5s,s'. (C.15) 



some useful definitions 



S'^{p;mj) = Y,{q{-p)i's{p))'{<'~'PL + w'^-^PR) = -S^{p;mf), (C.16) 

L« _ 

SL{jp;mf) = Y.^q{-p)i^s{p))'{<'~'PR + <-^PL) = S.^{p-mf). (C.17) 



We can, however, see that the one-loop heavy-light vertex correction (3.18) with the above 
Feynman rules is exactly same as Eq. (B.16): 

8T',,^,{p;m,) (C.18) 
^^U^)W^.^^(0, k)S±h{k)T Y,{<li-P - k)i^s{P + k))X^{p + k,p)s,sK{p;mj) 



x{-i)gT^\ - {S'^ip + k;mf) - S'^Xp + k;0)womf) 



cos ■ 



{k + 2p), 



{SwiP + k;mf) - S'^{p + k;0)womf) isin 



{k + 2p)^ 



0{p'^ ,pmf,m'j) 



8'^h±q{p; rriq). 



The expressions of the light quark wave function renormalization and the light-light vertex 
correction are also not altered by the change of the action convention. 
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D Power divergence structure in equations of motion 



In this appendix, we look into the power divergence structure in the equations of motion 
diagrammatically as the way in Sec. 4.4. In the analysis, we consider the heavy-light transi- 
tion amplitude of operators h±T^q for the light quark equation of motion and ^±r(7o^o)<J' 
for the static heavy quark equation of motion. Although these operators are identically 
vanished in the on-shell condition, each Feynman diagrams can have the divergences which 
are eventually canceled out. 

Equation of motion for light quark 

The light quark mass in the equation of motion is omitted, since this part does not cause 
the power divergence. Owing to the anisotropy of the static heavy-light system, the power 
divergent structure is different between temporal and spatial direction. The structure in 
the temporal direction is: 




- [0 + 0{j>a)] , 



[0 + Oig\pa)] , 



(D.l) 



(D.2) 



h±T{"/aDo)q 




'i±r(7oDo)<? 



-[0 + O(/,pa)], 



(D.3) 



'^^'V , = [+G6M + 0{9\pa)] , 



a 














p 





[0 + O{g\pa)] , 



/i±r(7oZ3o)9 

while in the spatial direction i (= 1,2,3) is: 



(D.4) 



(D.5) 



p 



1 



9 = - [0 + Oipa)] , 

h±r{r,D,)q « 




h±T{r,D,,)q 



^G5M + 0{g\pa) 



(D.6) 



(D.7) 
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In the temporal direction, the divergence comes from the loop diagram (D.4), but from the 
diagram (D.7) in the spatial direction. These divergences are totally cancelled out in their 
summation. 



Equation of motion for static quark 

Because the static quark action employed in this work does not use the symmetric co- 
variant derivative, we consider forward and backward derivatives separately to analyze the 
divergence structure. One definition of the covariant derivative, which is read directly from 
the static quark action (4.2), is 

h+l5o = h+{x + b)Ul{x) - h+{x), (D.ll) 
h^l5o = h^{x) -h^{x- d)Uo{x - 0). (D.12) 

This definition has following power divergence structure: 



''-±r(7o-Do)ij 



1 



- [GSM + 0(/)] 



(D.13) 




h±T{-foDo)q 



/i±r(7oDo)<7 



-[0 + O(/)], (D.14) 




p' fc3 k+ p' 



/i±r(7oDo)g 



-G6M + G^CpTi + 0(/) 



(D.15) 




/i±r(7oZ)o)g 



-G|-Ci.r4 + o(/) 



(D.16) 
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The tree diagram (D.13) generates a 1/a power divergence when the one-loop renormahzed 
on-shell condition for the static quark is imposed. The loop diagrams also contain power 
divergences. An interesting fact is that another 1/a power divergence T4, which is defined 
by Eq. (B.9), emerges in the diagrams (D.15) and (D.16). In totaL all 1/a divergences in 
the operator cancel, as is consistent with the equation of motion hD^ = 0. 
Another definition of the covariant derivative is 

h+l5'Q = h+{x) -h+{x-6)Uoix-6), (D.17) 
h-l5o = h_{x + 6)U^{x) - h_{x). (D.18) 

In this definition, we have the divergence structure 



hi Si □ M q = - [G5M + 0(/)] , (D.19) 

/i±r(7oBi)? « 




Again, the T4 divergences appear in the diagrams (D.21) and (D.22). Note that if we 
separately combine diagram (D.15) with (D.21) and (D.16) with (D.22) respectively, the 
divergences T4 vanish. The symmetric derivative does not induce divergences proportional 
to T4, which is consistent with the V ■ T symmetry. 

E Operator mixing in the static heavy and Wilson hght quark system 

In the paper we have used a domain-wall fermion action with chiral symmetry for the light 
quarks. In this appendix we briefly describe the pattern of operator mixing for Wilson 
fermions, for which there is explicit chiral symmetry breaking. 

Operator mixing for Wilson fermions with Wilson parameter r can conveniently be 
investigated by extending r to become a matrix R and by requiring consistency with the 
extended chiral symmetry, SU{Nf) x SU{Nf)- The chiral transformation of R matrix is 
deflned as 

R — > UlRUI, R^ UbR^uI, (E.l) 
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similar to Eq. (2.14). 

Analogous to Eq. (2.76), the matching of bilinear operators is 

(f|<^±r|i)cHQET = (f|^r^±r|i)LHQET, (E.2) 

with 















1 J±v = 


aJ±TD 




ZrM{r^)G + Zi^j^ir^y 




aJ±TM 



where each Z*{r'^) denotes functions of r^. For the matching of the AB = 2 four quark 



operators, we need following operators: 




Ol 




(E.4) 


On 


= 2[h-f;^q][h-fl:q]+A[hPRq][hPLq], 


(E.5) 


On 


= 2[7^7^g][7^7^<?] - A[hPRq][hPLq], 


(E.6) 


Or 


= [hll!q][hll!q], 


(E.7) 


Os 


= [hPLq][hPLq], 


(E.8) 


Ot 


= [hPRq][hPRq], 


(E.9) 


Old 


= [7^7^(7 •B)(/][7^7^g], 


(E.IO) 


OsD 


= [hPLh-D)q][hPLq], 


(E.ll) 


Ond 


= 2[7^7^(7 • ^)q][h^f:q] + 4[7^Pr(7 • ^)q][hPLq], 


(E.12) 




= 2[7^7^(7 • t)q][h^f:q] - A[hPR{f • t)q][hPLq], 


(E.13) 


Ond' 


= 2[7^7^(7 • t)q][hjl}q] + 4[7^Pi(7 • D)q][hPRq], 


(E.14) 


Ond' 


= 2[7^7^(7 • t)q][hjl}q] - 4[hPLh ■ D)q][hPRq], 


(E.15) 


Ord 


= [hj^{j-t)q][hj^q], 


(E.16) 


Otd 


= [hPR{f^)q][hPRq], 


(E.17) 


Olm 


= mq\hl^qWl^q]: 


(E.18) 


OsM 


= mg[hPLq]\hPLq], 


(E.19) 


Onm 


= {2[hj;!q][h^f:q]+4[hPRq][hPLq]) , 


(E.20) 


Onm 


= mg {2[hjl!q][h7f:q] - A[hPRq][hPLq]) , 


(E.21) 


Orm 


= mq[h-fl^q][h-fj^q], 


(E.22) 


Otm 


= mg\hPRq][hPRq], 


(E.23) 


and the matching relations are modified from Eqs. (2.84) and (2.85) to 






(f Ol i)cHQET = (f ^^fOl i)LHQET, 


(E.24) 




(f|0L + 40s|i)cHQET = (f|^2''02|i)LHQET, 


(E.25) 
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with 



)r 

ZLD{r^)r 
ZND{r^) 
ZND'{r^)r^ 
ZRD{r^)r 
ZLM{r^)r 
ZNM{r'') 
ZRM{r^)r 

ZL+As{r^) 

ZR+iT{r^y 
ZLD+ASD{r'^)r 
Zjmir^) 

ZRD+iTD{r'^)r 
ZLM+ASM{r'^)r 

ZRM+iTM{r'^)r 



Ol 
On 
Or 
aOLD 
aOND 
aO^D' 
aORD 
aOiM 
o-Onm 
iOrm 



(E.26) 



Ol + 4O5 



On 



Or + 40t 

a{OLD + 40sD) 



''Ond' 



a{ORD + 40™) 
a{OLM + 40sm) 

^OnJT 
a{ORM + 4:0tm) 



(E.27) 



The one-loop matching factor of the static heavy and clover- Wilson light quark system 
was calculated by Ishikawa et al. [8], in which the 0{pa) correction was included but the 
0{ma) correction was not. Note that operators Ond' and Ord in Eq. (E.26) are not seen 
in Ref. [8], because these operators arise at 0{g^). We can partly see that these operators 
actually arise from cross terms of 0{g'^) contributions in Eqs. (E.2) and (E.3), which leads 
to 0{g^). 

In the domain-wall quark formalism, similar intrinsic chiral symmetry breaking occurs 
for finite extent of the fifth dimension Ls- The size of the breaking effect is characterized 
by the residual mass ciTTi^Qgj which is suppressed as Lg increases. In principle, similar 
analysis for the operator mixing for the domain- wall fermion with finite Lg would be possible 
by introducing extended chiral symmetry for anires 7^ similar to Eq. (E.l). In actual 
numerical simulations using the domain-wall fermion, the residual mass is well controlled, 
arrires ~ O(10~^), which is negligible in practice. 
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F Tables of numerical values of domain-wall light quark part 

The values in this appendix are all cited from Refs. [32, 55]. 



1^5 


Plaquette 


Symanzik 


Iwasaki 


DBW2 


0.1 


-51 


048195 


-40 


074665 


-25 


740748 


-11 


263725 


0.2 


-50 


744997 


-39 


788806 


-25 


489672 


-11 


080847 


0.3 


-50 


488509 


-39 


548978 


-25 


282889 


-10 


936985 


0.4 


-50 


266419 


-39 


342775 


-25 


107854 


-10 


819677 


0.5 


-50 


072616 


-39 


163976 


-24 


958176 


-10 


722539 


0.6 


-49 


903778 


-39 


009135 


-24 


830209 


-10 


641842 


0.7 


-49 


758186 


-38 


876383 


-24 


721839 


-10 


575304 


0.8 


—49 


635208 


—38 


764914 


—24 


631970 


— 10 


521558 


0.9 


-49 


535083 


-38 


674758 


-24 


560281 


-10 


479897 


1.0 


-49 


458848 


-38 


606702 


-24 


507131 


-10 


450158 


1.1 


-49 


408369 


-38 


562312 


-24 


473568 


-10 


432689 


1.2 


-49 


386472 


-38 


544048 


-24 


461408 


-10 


428393 


1.3 


-49 


397179 


-38 


555483 


-24 


473426 


-10 


438846 


1.4 


-49 


446110 


-38 


601678 


-24 


513680 


-10 


466525 


1.5 


-49 


541138 


-38 


689808 


-24 


588063 


-10 


515219 


1.6 


-49 


693506 


-38 


830222 


-24 


705264 


-10 


590776 


1.7 


-49 


919841 


-39 


038393 


-24 


878568 


-10 


702549 


1.8 


-50 


246184 


-39 


338831 


-25 


129527 


-10 


866515 


1.9 


-50 


717592 


-39 


774527 


-25 


497027 


-11 


113350 



Table 4. Numerical values of Suj. The values are cited fronr Ref. [32]. 
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Ms 


riaquette 


oymanziK 


J 

Iwasaki 






0.1 


13.1613(11) 


9.5458(11) 


4.6519(11) 


-0.5152(11) 




0.2 


13.0109(11) 


9.4009(11) 


4.5193(11) 


-0.6190(11) 




0.3 


10 0000/1 1\ 

12.8838(11) 


9.2791(11) 


4.4093(11) 


-0.7020(11) 




0.4 


12.7739(11) 


9.1743(11) 


4.3158(11) 


-0.7704(11) 




0.5 


12.6781(11) 


9.0834(11) 


4.2354(11) 


-0.8276(11) 




0.6 


12.5948(11) 


9.0046(11) 


4.1665(11) - 


-0.8754(11) 




0.7 


12.5230(11) 


8.9370(11) 


4.1079(11) 


-0.9151(11) 




0.8 


12.4625(11) 


8.8803(11) 


4.0593(11) 


-0.9474(11) 




0.9 


12.4133(11) 


8.8344(11) 


4.0204(11) 


-0.9725(11) 




1.0 


12.3760(11) 


8.7998(11) 


3.9915(11) 


-0.9905(11) 




1.1 


12.3513(11) 


8.7773(11) 


3.9731(11) 


-1.0012(11) 




1.2 


12.3408(11) 


8.7681(11) 


3.9664(11) 


-1.0041(11) 




1.3 


12.3464(11) 


8.7740(11) 


3.9727(11) 


-0.9981(11) 




1.4 


1 Qvnc/' 1 1 \ 
iz.o ( Uo^ii j 


0. ( y ( o(ii j 


Q on /I Q/i 1 \ 

o.yy4o(iij - 


-U.yoio(^ii j 




1.5 


12.4179(11) 


8.8426(11) 


4.0343(11) - 


-0.9531(11) 




1.6 


12.4932(11) 


8.9141(11) 


4.0974(11) - 


-0.9085(11) 




1.7 


12.6050(11) 


9.0200(11) 


4.1905(11) - 


-0.8427(11) 




1.8 


12.7661(11) 


9.1726(11) 


4.3249(11) - 


-0.7469(11) 




1.9 


12.9989(11) 


9.3936(11) 


4.5209(11) - 


-0.6041(11) 




Ms 


Plaquette 


Symanzik 


V 

Iwasaki 


DBW2 


0.1 


-3.8182(35) 


-3.4782(35) 


-2.7691(35) 


-1.1915(35) 


0.2 


-4.6031(37) 


-4.2403(36) 


-3.4827(37) 


-1.7967(36) 


0.3 


-5.2388(34) 


-4.8520(35) 


-4.0444(34) 


-2.2502(34) 


0.4 


-5.7981(36) 


-5.3863(36) 


-4.5265(36) 


-2.6240(36) 


0.5 


-6.3138(36) 


-5.8754(36) 


-4.9616(36) 


-2.9496(36) 


0.6 


-6.8035(36) 


-6.3371(36) 


-5.3668(36) 


-3.2439(36) 


0.7 


-7.2795(35) 


-6.7834(36) 


-5.7538(35) 


-3.5176(35) 


0.8 


-7.7511(35) 


-7.2234(35) 


-6.1314(35) 


-3.7789(35) 


0.9 


-I 


^.2259(36) 


-7.6643(37) 


-6.5064(36) 


-4.0335(36) 


1.0 


-I 


5.7112(36) 


-8.1136(35) 


-6.8854(36) 


-4.2871(35) 


1.1 


-9.2169(19) 


-8.5801(19) 


-7.2772(19) 


-4.5470(19) 


1.2 


-9.7463(30) 


-9.0671(30) 


-7.6839(30) 


-4.8141(30) 


1.3 


-10.3140(34) 


-9.5888(33) 


-8.1186(34) 


-5.0994(34) 


1.4 


-10.9324(33) 


-10.1564(37) 


-8.5918(33) 


-5.4109(34) 


1.5 


-11.6187(37) 


-10.7872(37) 


-9.1182(37) 


-5.7606(37) 


1.6 


-12.3986(37) 


-11.5052(34) 


-9.7205(37) 


-6.1668(36) 


1.7 


-13.3112(37) 


-12.3479(37) 


-10.4338(37) 


-6.6591(37) 


1.8 


-14.4234(37) 


-13.3811(36) 


-11.3197(37) 


-7.2913(38) 


1.9 


-15.8762(35) 


-14.7426(37) 


-12.5115(35) 


-8.1852(35) 



Table 5. Numerical values of / and v. The values are produced using the data in Ref. [32]. 
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M5 




M5 


df 


0.05 


—0.028209(17) 


1.05 


—0.0037322(49) 


0.10 


—0.023030(14) 


1.10 


—0.0032266(45) 


0.15 


—0.020075(13) 


1.15 


—0.0027422(43) 


0.20 


—0.017981(11) 


1.20 


—0.0022896(41) 


0.25 


— O.Ulo332(ll) 


i.ZO 


r\ r\ni cvi o/on\ 

—0.0018718(39) 


O.oO 


—0.0149658(97) 


1.30 


—0.0014846(37) 


not; 

0.35 


-0.0137846(96) 


1.35 


—0.0011441(36) 


0.40 


-0.0127360(87) 


1.40 


-0.0008650(35) 


0.45 


—0.0117805(83) 


1.45 


—0.0006560(34) 


0.50 


—0.0109029(79) 


1.50 


—0.0005360(33) 


0.55 


—0.0100864(76) 


1.55 


—0.0005236(33) 


0.60 


—0.0093153(73) 


1.60 


—0.0006566(34) 


0.65 


—0.0085901(67) 


1.65 


—0.0009842(36) 


0.70 


—0.0078955(64) 


1.70 


r\ r\r\'\ f '~7r\ a / '~tf\\ 

—0.0015704(39) 


7^ 


—0 0n799q4('63'l 


1 75 


—0 n09590n('43'i 


0.80 


-0.0065893(62) 


1.80 


-0.0040140(51) 


0.85 


-0.0059729(55) 


1.85 


-0.0063612(63) 


0.90 


-0.0053788(55) 


1.90 


-0.0102019(82) 


0.95 


-0.0048118(51) 


1.95 


-0.017386(12) 


1.00 


-0.0042608(51) 







Table 6. Numerical values oi df. The values are cited from Ref. [55]. 
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G Tables of numerical values of integrals 



To perform the numerical integration, we use the Monte Carlo integration routine VEGAS 
[62]. 



gluon action 


unsmear 


7^^l' 

APE 


-PL 

HYPl 


HYP2 


Plaquette 


0.0 


-8.155358(22) 


-7.768507(29) 


-12.346005(55) 


Symanzik 


-0.595480(4) 


-8.157309(24) 


-8.061365(30) 


-11.862191(52) 


Iwasaki 


-2.620439(12) 


-8.610915(27) 


-8.769236(31) 


-11.574859(51) 


DBW2 


-7.020924(38) 


-10.287034(43) 


-10.528146(45) 


-12.174154(55) 



Table 7. Numerical values of integrals TZ"^^ 



gluon action 


unsmear 


^3 

APE 


-PL 

HYPl 


HYP2 


Plaquette 


0.0 


-0.166666(1) 


-0.179744(1) 


-0.199563(1) 


Symanzik 


-0.033703(0) 


-0.172377(0) 


-0.183290(1) 


-0.201904(1) 


Iwasaki 


-0.088361(0) 


-0.182888(1) 


-0.190549(1) 


-0.206241(1) 


DBW2 


-0.158892(1) 


-0.201084(1) 


-0.204909(1) 


-0.214596(1) 



Table 8. Numerical values of integrals Tg^ 



gluon action 




6M 




unsmear 


APE 


HYPl 


HYP2 


Plaquette 


19.95473(10) 


6.79531(11) 


5.76272(11) 


4.19788(13) 


Symanzik 


17.29363(10) 


6.34435(11) 


5.48274(11) 


4.01300(13) 


Iwasaki 


12.97799(11) 


5.51444(11) 


4.90959(11) 


3.67056(12) 


DBW2 


7.40912(12) 


4.07778(12) 


3.77575(12) 


3.01091(13) 



Table 9. Numerical values of SM. 



gluon action 




APE 


e 




unsmear 


HYPl 


HYP2 


Plaquette 


24.48049(10) 


3.16571(11) 


2.51997(12) 


-3.62237(14) 


Symanzik 


21.22391(10) 


2.71280(12) 


1.94714(12) 


-3.32344(14) 


Iwasaki 


14.88332(11) 


1.42929(12) 


0.66611(12) 


-3.37854(14) 


DBW2 


4.91396(13) 


-1.68349(13) 


-2.22663(13) 


-4.63749(14) 



Table 10. Numerical values of e. 
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1\ /T 

Ms 


7-PL-div 


tPL 


o-PL-div 




0.1 


—1.317357(54) 


—4.034531(33) 


—2.55733(63) 


13.49853(12) 


0.2 


— 1.490882(49) 


—3.872305(25) 


—0.47980(56) 


1 1 0'70000/'00\ 

11.072832(82) 


0.3 


—1.648560(46) 


—3.738575(21) 


0.73748(38) 


9.680365(61) 


0.4 


1 '1Q n r\r\r\ ( a a\ 

— 1.786090(44) 


/^oooi /i/ioN 

—3.620014(18) 


1.62142(33) 


8.689872(48) 


0.5 


— 1.915172(44) 


—3.510747(16) 


n 00001 /oo\ 

2.33301(29) 


'1 00*70 /I 1 / 00^ 

7.907941(39) 


O.b 


—2.039462(43) 


A r\''7 A A r\ ^ '\ a\ 

—3.407449(14) 


2.94407(27) 


7.249261(31) 


0.7 


-2.161204(43) 


00'7000/'10\ 

—3.307920(12) 


/iooor/oo\ 

3.49325(30) 


n C)C\'\ 0/01 \ 

6.668012(21) 


0.8 


—2.282065(43) 


OIOCTOO/ION 

—3.210533(10) 


/I 00 /I 1 /oo\ 

4.00481(32) 


6.135881(16) 


0.9 


—2.403557(45) 


—3.113948(9) 


4.49564(27) 


5.633217(11) 


1.0 


—2.527097(45) 


—3.016961(6) 


4.97879(29) 


5.144634(4) 


1.1 


—2.654088(46) 


—2.918368(8) 


5.46667(30) 


4.656455(10) 


1.2 


—2.785922(47) 


—2.816917(9) 


r" 0^1 1 yi/oo\ 

5.97114(39) 


4.154848(13) 


1.3 


—2.924429(48) 


— 2.711123(10) 


6.50702(41) 


3.623851(16) 


1 4 


—3 071 769('4Q'l 


— 9 59Q1 fion 1 


7 0Q1 Q^Md^l 


3 049fifi3('99') 


1.5 


-3.230735(51) 


-2.478608(12) 


7.75540(40) 


2.380852(30) 


1.6 


-3.405303(53) 


-2.346024(14) 


8.53794(43) 


1.588916(37) 


1.7 


-3.601408(58) 


-2.196048(16) 


9.52090(47) 


0.574365(48) 


1.8 


-3.829124(59) 


-2.019452(20) 


10.88105(65) 


-0.871646(65) 


1.9 


-4.109490(66) 


-1.796311(27) 


13.18814(87) 


-3.42206(10) 



Ms 


^PL-div 
x' 


^PL-div 


0.1 


-0.03163(56) 


-41.6877(16) 


0.2 


1.23722(41) 


-19.57253(70) 


0.3 


1.99022(35) 


-11.74501(46) 


0.4 


2.54399(30) 


-7.55285(36) 


0.5 


2.99506(27) 


-4.82466(43) 


0.6 


3.38651(26) 


-2.82047(29) 


0.7 


3.74227(26) 


-1.21451(20) 


0.8 


4.07642(27) 


0.16622(17) 


0.9 


4.39987(29) 


1.42852(17) 


1.0 


4.72007(34) 


2.64870(20) 


1.1 


5.04472(32) 


3.89320(23) 


1.2 


5.38204(35) 


5.23129(33) 


1.3 


5.74066(34) 


6.75210(44) 


1.4 


6.13289(37) 


8.58591(35) 


1.5 


6.57496(39) 


10.95370(62) 


1.6 


7.09382(45) 


14.28337(50) 


1.7 


7.73863(46) 


19.55246(79) 


1.8 


8.61765(53) 


29.65130(89) 


1.9 


10.07462(69) 


58.8508(17) 



Table 11. Numerical values of integrals Z^L-div^ jPL^ J^^"^^^ , Jl^ , 3^}^"^" and /C^ 
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unsmear 


7-5! -PL 
APE HYPl 




0.1 


0.0 


0.407534(12) 


0.393862(14) 


0.621219(25) 


(J. 2 


0.0 


0.418636(13) 


0.405266(15) 


0.639854(26) 


(J. 3 


0.0 


0.430385(13) 


0.417365(15) 


0.659656(27) 


(J. 4 


0.0 


0.442837(13) 


0.430221(16) 


0.680726(28) 


(J. 5 


0.0 


0.456053(14) 


0.443901(16) 


0.703179(29) 


(J.D 


0.0 


0.470102(14) 


0.458481(17) 


0.727146(30) 


(J. 7 


0.0 


0.485064(15) 


0.474051(17) 


0.752777(31) 


0.8 


0.0 


0.501027(16) 


0.490709(18) 


0.780244(32) 


0.9 


0.0 


0.518093(16) 


0.508571(19) 


0.809745(34) 


1.0 


0.0 


0.536381(17) 


0.527771(20) 


0.841510(35) 


1.1 


0.0 


0.556028(18) 


0.548465(21) 


0.875811(37) 


1.2 


0.0 


0.577192(18) 


0.570837(22) 


0.912966(39) 


1.3 


0.0 


0.600063(19) 


0.595104(23) 


0.953355(41) 


1 A 


U.U 


0.624863(20) 


0.621527(24) 


u.yy 1 4oo(^4o ) 


1.5 


0.0 


0.651860(22) 


0.650423(26) 


1.045763(46) 


1.6 


0.0 


0.681381(23) 


0.682178(27) 


1.099024(49) 


1.7 


0.0 


0.713828(25) 


0.717277(30) 


1.158079(53) 


1.8 


0.0 


0.749705(27) 


0.756331(32) 


1.224021(58) 


1.9 


0.0 


0.789656(30) 


0.800130(36) 


1.298275(65) 




Ms 


unsmear 


^x 

APE HYPl 


HYP2 


0.1 


0.0 


-2.037596(54) 


-2.365827(68) 


-4.13518(14) 


0.2 


0.0 


-2.100353(55) 


-2.440650(69) 


-4.26767(14) 


0.3 


0.0 


-2.163706(56) 


-2.515804(70) 


-4.40038(14) 


0.4 


0.0 


-2.228348(57) 


-2.592118(80) 


-4.53494(14) 


0.5 


0.0 


-2.294787(58) 


-2.670352(81) 


-4.67253(15) 


0.6 


0.0 


-2.363474(59) 


-2.751008(83) 


-4.81415(15) 


0.7 


0.0 


-2.434800(67) 


-2.834629(84) 


-4.96078(15) 


0.8 


0.0 


-2.509381(61) 


-2.921776(85) 


-5.11343(16) 


0.9 


0.0 


-2.587563(62) 


-3.013065(87) 


-5.27319(16) 


1.0 


0.0 


-2.669967(63) 


-3.109187(89) 


-5.44131(16) 


1.1 


0.0 


-2.757255(64) 


-3.210948(90) 


-5.61923(16) 


1.2 


0.0 


-2.850204(65) 


-3.319302(92) 


-5.80866(17) 


1.3 


0.0 


-2.949756(66) 


-3.435408(93) 


-6.01168(17) 


1.4 


0.0 


-3.057067(68) 


-3.560702(95) 


-6.23089(17) 


1.5 


0.0 


-3.173600(69) 


-3.697012(97) 


-6.46961(18) 


1.6 


0.0 


-3.301186(79) 


-3.84673(10) 


-6.73217(18) 


1.7 


0.0 


-3.442518(81) 


-4.01310(10) 


-7.02452(19) 


1.8 


0.0 


-3.601161(83) 


-4.20078(10) 


-7.35520(19) 


1.9 


0.0 


-3.782775(87) 


-4.41709(11) 


-7.73773(20) 


Table 12. Numerical values of integ 


rals ^x''^^ and J^'^^^ (Plaquette) 
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unsmear 


APE 


— PT, 

HYPl 


HYP 2 


0.1 


0.0 


-4.283384(10) 


-3.984201(16) 


-6.147964(32) 


0.2 


0.0 


-4.209373(9) 


-3.895649(14) 


-5.990706(28) 


0.3 


0.0 


-4.137273(8) 


-3.810173(13) 


-5.839628(28) 


0.4 


0.0 


-4.065997(7) 


-3.726307(12) 


-5.692018(25) 


0.5 


0.0 


-3.994823(7) 


-3.643141(12) 


-5.546171(24) 


0.6 


0.0 


-3.923205(6) 


-3.559974(11) 


-5.400810(23) 


0.7 


0.0 


-3.850673(6) 


-3.476212(11) 


-5.254859(23) 


0.8 


0.0 


-3.776781(6) 


-3.391307(11) 


-5.107324(22) 


0.9 


0.0 


-3.701087(7) 


-3.304709(11) 


-4.957219(22) 


1.0 


0.0 


-3.623098(7) 


-3.215843(11) 


-4.803507(22) 


1.1 


0.0 


-3.542290(6) 


-3.124055(12) 


-4.645032(22) 


1.2 


0.0 


-3.458053(7) 


-3.028608(12) 


-4.480485(23) 


1.3 


0.0 


-3.369660(7) 


-2.928625(12) 


-4.308282(23) 


1.4 


0.0 


-3.276227(7) 


-2.823037(11) 


—4.126494(24) 


1.5 


0.0 


-3.176639(8) 


-2.710462(12) 


-3.932684(23) 


1.6 


0.0 


-3.069451(8) 


-2.589120(12) 


-3.723605(24) 


1.7 


0.0 


-2.952709(9) 


-2.456580(13) 


-3.494856(24) 


1.8 


0.0 


-2.823657(13) 


-2.309344(14) 


-3.240068(26) 


1.9 


0.0 


-2.677956(12) 


-2.141922(19) 


-2.949198(38) 










x' 

HYPl 




Ms 


unsmear 


APE 


HYP2 


0.1 


0.0 


-2.037576(16) 


-2.224503(20) 


-3.773845(37) 


0.2 


0.0 


-2.100332(16) 


-2.295041(21) 


-3.895183(38) 


0.3 


0.0 


-2.163685(16) 


-2.365678(21) 


-4.016136(38) 


0.4 


0.0 


-2.228325(17) 


-2.437291(21) 


-4.138302(39) 


0.5 


0.0 


-2.294764(17) 


-2.510505(22) 


-4.262808(40) 


0.6 


0.0 


-2.363451(17) 


-2.585858(22) 


-4.390612(40) 


0.7 


0.0 


-2.434827(18) 


-2.663872(22) 


-4.522630(41) 


0.8 


0.0 


-2.509355(18) 


-2.745082(23) 


-4.659803(42) 


0.9 


0.0 


-2.587537(18) 


-2.830075(23) 


-4.803152(42) 


1.0 


0.0 


-2.669943(18) 


-2.919507(24) 


-4.953826(43) 


1.1 


0.0 


-2.757232(19) 


-3.014145(24) 


-5.113162(44) 


1.2 


0.0 


-2.850182(19) 


-3.114897(24) 


-5.282749(44) 


1.3 


0.0 


-2.949748(26) 


-3.222887(33) 


-5.464530(45) 


1.4 


0.0 


-3.057059(26) 


-3.339448(34) 


-5.660927(46) 


1.5 


0.0 


-3.173594(27) 


-3.466349(34) 


-5.875094(63) 


1.6 


0.0 


-3.301250(28) 


-3.605891(35) 


-6.111047(64) 


1.7 


0.0 


-3.442586(28) 


-3.761213(36) 


-6.374430(49) 


1.8 


0.0 


-3.601226(22) 


-3.936824(28) 


-6.673536(51) 


1.9 


0.0 


-3.782849(23) 


-4.139930(29) 


-7.021409(54) 


Table 13. Numerical values of inte 


grals J.^^-^^ and J^r^^ (Plaquette) 
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Ms 


unsmear 


APE 


HYPl 


HYP2 


0.1 


0.0 


17.853525(31) 


15.988954(48) 


24.027654(94) 


0.2 


0.0 


7.811064(21) 


6.811607(32) 


10.035191(61) 


0.3 


0.0 


4.392120(18) 


3.679576(26) 


5.251770(48) 


0.4 


0.0 


2.621649(16) 


2.052067(22) 


2.760193(42) 


0.5 


0.0 


1.502881(15) 


1.019094(20) 


1.174087(38) 


0.6 


0.0 


0.701465(14) 


0.275394(19) 


0.028250(36) 


0.7 


0.0 


0.071833(15) 


-0.312087(19) 


-0.880199(34) 


0.8 


0.0 


-0.461561(15) 


-0.812625(19) 


-1.657020(47) 


0.9 


0.0 


-0.944408(16) 


-1.268104(20) 


-2.366378(37) 


1.0 


0.0 


-1.408989(17) 


-1.708448(21) 


-3.054217(38) 


1.1 


0.0 


-1.882758(18) 


-2.159268(22) 


-3.760153(41) 


1.2 


0.0 


-2.394296(18) 


-2.647479(24) 


-4.526036(43) 


1.3 


0.0 


-2.979644(20) 


-3.207259(25) 


-5.405258(46) 


1.4 


0.0 


-3.692200(28) 


-3.889483(37) 


-6.477516(67) 


1.5 


0.0 


-4.622672(30) 


-4.780726(40) 


-7.878588(74) 


1.6 


0.0 


-5.947873(34) 


-6.049948(45) 


-9.873656(83) 


1.7 


0.0 


-8.075885(38) 


-8.087229(52) 


-13.075106(70) 


1.8 


0.0 


-12.226002(47) 


-12.058228(48) 


-19.312790(89) 


1.9 


0.0 


-24.487147(52) 


-23.783853(74) 


-37.72516(14) 




Table 14. 


Numerical values of integrals JC^ 


(Plaquctte). 



- 52 - 





unsmear 




-PL 

M Yr 1 




0.1 


0.102249(1) 


0.420466(9] 


A 'X r\r\c\'\ ( 'X "X \ 

0.410991(11] 


r'OOyloo/io\ 

0.590402(19) 


0.2 


0.105228(1] 


/100100/'10\ 

0.432189(10) 


/100000/ll\ 

0.422982(11] 


^001 o'7/oo\ 

0.608137(20] 


L).3 


0.108389(1] 


/1/I/1/^00/10\ 

0.444603(10] 


/ior'70o/io\ 

0.435703(12) 


^0/^0'70/00\ 

0.626973(20] 


(J. 4 


0.111749(2] 


0.457769(10) 


/1/1001/^/lO\ 

0.449216(12) 


A'1C\r\ A ( ^'X\ 

0.647004(21] 


U.o 


0.115o2d(2] 


/1'71'V/10/'ll\ 

0.471749(11) 


/i^ocroo/io\ 

0.463592(12) 


^/^000/1/00\ 

0.668334(22] 


U.D 


0.119140(2] 


/10/^/^10/'ll\ 

0.486618(11) 


/1'70000/10\ 

0.478909(13) 


^oioor/oo\ 

0.691085(23] 


L).7 


r\ 10001 /i/o\ 

0.123214(2] 


r'oo/ir'o/ii\ 

0.502459(11) 


/ioror'o/io\ 

0.495258(13) 


'71 CTOO/l/O/lX 

0.715394(24] 


0.8 


io^p'^o/o\ 

0.127578(2] 


0.519371(12) 


V 'X A A ( "X A\ 

0.512744(14) 


'1 A "X /101 ( C\T'\ 

0.741421(25] 


0.9 


iooor'o/o\ 

0.132258(2] 


0.537456(12) 


r'oi /I00/1 a\ 

0.531482(14) 


0.769343(26] 


1.0 


1 o'^ooo/oN 

0.137290(2] 


0.556844(13) 


0.551610(15) 


'^ooo'^r'/o'^N 

0.799375(27] 


1.1 


1 /io'^i'^/oN 

0.142717(2] 


v '1*1 n c^'X /1 o\ 

0.577681(13) 


r''^oooo/i/^\ 

0.573292(16) 


ooi^/^o/ooN 

0.831763(28] 


1.2 


1 /lOCTOO/oX 

0.148589(2] 


/^ooioo/i a\ 

0.600138(14) 


0.596714(17) 


o/^/^ooi /"ooX 

0.866801(29] 


1.3 


1 cr/io/^o/o\ 

0.154962(2] 


/^O/iyii'7/1 r\ 

0.624417(15) 


/^00100/l'7\ 

0.622103(17) 


00/100'7/'01\ 

0.904837(31] 


1 4 


1fi1Q07('9'l 


fi5075Q('1 6") 


fidQT^sn 8") 


n 94698Ql"^3'l 


1.5 


0.169510(2) 


0.679455(17) 


0.679915(19) 


0.991668(35) 


1.6 


0.177874(3) 


0.710859(18) 


0.713064(21) 


1.041601(37) 


1.7 


0.187131(3) 


0.745409(19) 


0.749676(22) 


1.096875(40) 


1.8 


0.197446(3) 


0.783657(20) 


0.790383(24) 


1.158490(44) 


1.9 


0.209036(3) 


0.826311(23) 


0.836005(27) 


1.227747(48) 



Ms 


unsmear 


APE 


-PL 

HYPl 


HYP2 


0.1 


-0.553407(3) 


-2.390646(52) 


-2.708685(65) 


-4.34123(13) 


0.2 


-0.571158(3) 


-2.465179(53) 


-2.794633(67) 


-4.48021(14) 


0.3 


-0.589055(3) 


-2.540081(54) 


-2.880630(68) 


-4.61890(12) 


0.4 


-0.607304(3) 


-2.616213(55) 


-2.967729(69) 


-4.75889(13) 


0.5 


-0.626058(3) 


-2.694200(56) 


-3.056681(70) 


-4.90152(13) 


0.6 


-0.645451(3) 


-2.774581(57) 


-3.148126(72) 


-5.04784(13) 


0.7 


-0.665616(3) 


-2.857880(58) 


-3.242680(73) 


-5.19885(13) 


0.8 


-0.686693(3) 


-2.944636(59) 


-3.340975(74) 


-5.35558(13) 


0.9 


-0.708836(3) 


-3.035439(60) 


-3.443695(75) 


-5.51912(14) 


1.0 


-0.732219(3) 


-3.130949(61) 


-3.551606(76) 


-5.69062(16) 


1.1 


-0.757046(3) 


-3.231926(62) 


-3.665600(77) 


-5.87172(16) 


1.2 


-0.783559(3) 


-3.339278(63) 


-3.786732(79) 


-6.06401(16) 


1.3 


-0.812055(4) 


-3.454095(64) 


-3.916224(90) 


-6.26956(16) 


1.4 


-0.842901(4) 


-3.577726(65) 


-4.055774(92) 


-6.49095(17) 


1.5 


-0.876565(4) 


-3.711877(66) 


-4.207342(94) 


-6.73145(17) 


1.6 


-0.913662(4) 


-3.858768(68) 


-4.373567(96) 


-6.99535(17) 


1.7 


-0.955035(4) 


-4.021407(69) 


-4.558045(98) 


-7.28853(18) 


1.8 


-1.001875(4) 


-4.203997(80) 


-4.76591(10) 


-7.61944(18) 


1.9 


-1.056084(4) 


-4.413346(83) 


-5.00527(10) 


-8.00142(19) 



Table 15. Numerical values of integrals and J!j^ (Symanzik). 
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-PL 




Ms 


unsmear 


MYPl 




0.1 


—0.950580(2) 


—4.515816(10) 


—4.306330(15) 


—6.168544(28) 


0.2 


—0.930449(1) 


—4.427096(9) 


—4.203564(13) 


—6.001124(25) 


0.3 


—0.910991(1) 


—4.341257(8) 


—4.104837(12) 


—5.841055(23) 


0.4 


—0.891870(1) 


—4.256863(7) 


—4.008368(11) 


—5.685305(21) 


0.5 


— 0.872877(1) 


—4.173009(6) 


—3.913053(10) 


—5.531974(20) 


O.D 


—0.853858(1) 


—4.089013(6) 


—3.818055(10) 


—5.379683(19) 


0.7 


—0.834676(1) 


—4.004291(5) 


—3.722670(9) 


—5.227248(19) 


0.8 


—0.815203(1) 


—3.918296(5) 


—3.626248(9) 


—5.073594(18) 


0.9 


—0.795315(1) 


— 3.830492(6) 


—3.528146(9) 


— 4.917661(18) 


1.0 


—0.774875(1) 


—3.740285(6) 


—3.427696(9) 


—4.758349(17) 


1.1 


— 0.753736(1) 


— 3.647050(6) 


— 3.324136(9) 


— 4.594419(17) 


1.2 


—0.731726(1) 


—3.550047(6) 


—3.216611(10) 


A A A A O A / 1 f\\ 

—4.424484(19) 


1.3 


—0.708640(1) 


— 3.448411(6) 


— 3.104112(10) 


—4.246869(20) 


1 A 
1.4 




— o.o4iUDy(^D ) 


— Z.yooooZ[l\J ) 


— 4.Uoyo24(^2U ) 


1.5 


-0.658175(1) 


-3.226676(7) 


-2.858835(10) 


-3.859845(21) 


1.6 


-0.630065(1) 


-3.103476(8) 


-2.722385(10) 


-3.644406(20) 


1.7 


-0.599338(1) 


-2.969070(9) 


-2.573161(12) 


-3.408438(21) 


1.8 


-0.565185(1) 


-2.820029(10) 


-2.407055(14) 


-3.145095(23) 


1.9 


-0.526356(1) 


-2.651005(10) 


-2.217538(15) 


-2.843478(25) 




Ms 


unsmear 


x' 

APE 


-PL 

HYPl 


HYP2 


0.1 


-0.518474(4) 


-2.321568(16) 


-2.516149(20) 


-3.934517(36) 


0.2 


-0.535073(4) 


-2.393807(16) 


-2.595982(20) 


-4.060531(36) 


0.3 


-0.551753(4) 


-2.466300(17) 


-2.675549(21) 


-4.185501(37) 


0.4 


-0.568716(4) 


-2.539902(17) 


-2.755889(21) 


-4.311167(38) 


0.5 


-0.586110(4) 


-2.615232(17) 


-2.837732(21) 


-4.438731(38) 


0.6 


-0.604063(4) 


-2.692823(17) 


-2.921698(22) 


-4.569195(39) 


0.7 


-0.622701(4) 


-2.773186(18) 


-3.008372(22) 


-4.703504(40) 


0.8 


-0.642157(5) 


-2.856852(18) 


-3.098354(22) 


-4.842612(40) 


0.9 


-0.662575(5) 


-2.944397(18) 


-3.192292(23) 


-4.987546(41) 


1.0 


-0.684118(5) 


-3.036465(19) 


-3.290912(23) 


-5.139455(42) 


1.1 


-0.706978(5) 


-3.133803(19) 


-3.395053(24) 


-5.299665(42) 


1.2 


-0.731382(5) 


-3.237294(19) 


-3.505709(24) 


-5.469750(43) 


1.3 


-0.757609(5) 


-3.348003(19) 


-3.624090(24) 


-5.651628(44) 


1.4 


-0.786005(5) 


-3.467249(20) 


-3.751695(25) 


-5.847688(44) 


1.5 


-0.817022(6) 


-3.596708(20) 


-3.890449(25) 


-6.060999(45) 


1.6 


-0.851228(6) 


-3.738542(21) 


-4.042855(26) 


-6.295587(46) 


1.7 


-0.889431(6) 


-3.895707(21) 


-4.212360(26) 


-6.557014(47) 


1.8 


-0.932788(6) 


-4.072415(22) 


-4.403936(27) 


-6.853349(48) 


1.9 


-0.983140(7) 


-4.275259(23) 


-4.625453(28) 


-7.197456(51) 



Table 16. Numerical values of integrals and J^} (Symanzik). 
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unsmear 


APE 


HYPl 


HYP2 


0.1 


3.849344(10) 


18.488439(31) 


17.061540(45) 


23.805797(84) 


0.2 


1.649442(7) 


7.986934(21) 


7.196058(30) 


9.838490(56) 


0.3 


0.898962(5) 


4.407619(18) 


3.826504(24) 


5.060109(45) 


0.4 


0.509197(5) 


2.551179(21) 


2.073722(29) 


2.568738(40) 


0.5 


0.261965(4) 


1.375696(15) 


0.959767(20) 


0.980994(36) 


0.6 


0.084081(4) 


0.531777(14) 


0.156666(19) 


-0.167351(34) 


0.7 


-0.056353(4) 


-0.132810(14) 


-0.478629(18) 


-1.078794(33) 


0.8 


-0.175939(4) 


-0.697215(15) 


-1.020603(25) 


-1.858881(33) 


0.9 


-0.284736(4) 


-1.209231(16) 


-1.514365(20) 


-2.571793(35) 


1.0 


-0.389905(4) 


-1.702809(17) 


-1.992117(21) 


-3.263345(37) 


1.1 


-0.497591(4) 


-2.206883(18) 


-2.481506(22) 


-3.973186(40) 


1.2 


-0.614251(5) 


-2.751683(19) 


-3.011611(23) 


-4.743187(42) 


1.3 


-0.748088(5) 


-3.375423(20) 


-3.619401(25) 


-5.626782(45) 


1.4 


-0.911311(6) 


-4.134800(21) 


-4.359904(27) 


-6.703685(49) 


1.5 


-1.124723(6) 


-5.126245(30) 


-5.326878(39) 


-8.109874(71) 


1.6 


-1.428914(7) 


-6.537665(34) 


-6.703115(44) 


-10.110511(79) 


1.7 


-1.917625(8) 


-8.802843(39) 


-8.910722(50) 


-13.318052(92) 


1.8 


-2.871027(11) 


-13.217740(47) 


-13.210794(62) 


-19.562256(84) 


1.9 


-5.688473(19) 


-26.253637(71) 


-25.900451(71) 


-37.97904(13) 



Table 17. Numerical values of integrals IC^ (Symanzik). 
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unsmear 




-PL 

H Yrl 




0.1 


0.249872(3) 


0.446767(7] 


/1ylOylOO/'7\ 

0.442480(7] 


r'r'r'o'7yi/io\ 

0.555974(12] 


0.2 


0.257249(3] 


^iro/^oo/'7\ 

0.459632(7] 


0.455554(8] 


crTo'71 0/ 1 o\ 

0.572719(12] 




0.265073(3] 


^1 '70 O/^ / 

0.473262(7] 


0.469417(8] 


croo^io^i/io\ 

0.590484(13] 


0.4 


0.273383(4] 


0.487718(7] 


/10/110^l/0\ 

0.484134(8] 


/^ooocro/io\ 

0.609353(13] 


0.5 


n 000010 / /i\ 

0.282219(4] 


^000*71 { 'l\ 

0.503071(7] 


0.499776(8] 


/^oo^iit/i a\ 

0.629417(14] 


O.D 


r\ 001 ^00/ /I \ 

0.291629(4] 


crioooo/o\ 

0.519399(8] 


ri/^^iocr/o\ 

0.516425(9] 


/^cro'70o/i a\ 

0.650783(14] 


0.7 


0.301668(4] 


cro/^'70o/o\ 

0.536793(8] 


ro/ii'7/^/o\ 

0.534176(9] 


/^'70cr'7o/i r\ 

0.673573(15] 


0.8 


0.312397(4] 


0.555354(8] 


r'r'oior'/o\ 

0.553135(9] 


r'o^oor'/ir'\ 

0.697925(16] 


0.9 


000000 / /iN 

0.323888(4] 


0.575201(9] 


r'^0/io/^/"io\ 

0.573426(10] 


^0/iooi/ir'\ 

0.724001(16] 


1.0 


oo/^oor'//i\ 

0.336225(4] 


0.596472(9] 


0.595194(10] 


^r'l 000/1 ^\ 

0.751989(17] 


1.1 


0.349505(5] 


/^'iooor'/o\ 

0.619326(9] 


0.618606(11] 


^001 0*7/1 o\ 

0.782107(18] 


1.2 


o/^oo/io/cr\ 

0.363843(5] 


0.643951(10] 


/^/100/^0/l 1 \ 

0.643862(11] 


01 yi/^ir/io\ 

0.814615(18] 


1.3 


0.379375(5] 


0.670570(10] 


r"i'\ io^/io\ 

0.671196(12] 


0/10000/lO\ 

0.849823(19] 


1 4 


3Qfi9fi7('5'l 


n fiQQ44Q('1 1 


700891 i^ 9") 


n 888101 (9()\ 


1.5 


0.414717(6) 


0.730910(11) 


0.733292(13) 


0.929902(21) 


1.6 


0.434971(6) 


0.765348(12) 


0.768820(14) 


0.975784(23) 


1.7 


0.457337(6) 


0.803254(13) 


0.808001(15) 


1.026442(24) 


1.8 


0.482204(7) 


0.845247(14) 


0.851503(16) 


1.082764(26) 


1.9 


0.510080(8) 


0.892127(15) 


0.900194(17) 


1.145901(29) 



Ms 


unsmear 


APE 


-PL 

HYPl 


HYP2 


0.1 


-1.775035(10) 


-3.268508(50) 


-3.562292(62) 


-4.94892(11) 


0.2 


-1.831856(10) 


-3.371661(50) 


-3.675559(63) 


-5.10681(11) 


0.3 


-1.888377(10) 


-3.474103(51) 


-3.787699(64) 


-5.26262(12) 


0.4 


-1.945338(11) 


-3.577181(52) 


-3.900243(65) 


-5.41855(12) 


0.5 


-2.003257(11) 


-3.681825(53) 


-4.014239(66) 


-5.57608(12) 


0.6 


-2.062569(11) 


-3.788808(54) 


-4.130550(67) 


-5.73644(12) 


0.7 


-2.123681(11) 


-3.898850(55) 


-4.249976(68) 


-5.90075(12) 


0.8 


-2.187011(12) 


-4.012679(56) 


-4.373320(69) 


-6.07012(12) 


0.9 


-2.253006(12) 


-4.131072(56) 


-4.501438(70) 


-6.24573(13) 


1.0 


-2.322167(12) 


-4.254895(57) 


-4.635279(71) 


-6.42890(13) 


1.1 


-2.395072(12) 


-4.385139(58) 


-4.775936(73) 


-6.62112(13) 


1.2 


-2.472405(12) 


-4.522975(59) 


-4.924695(74) 


-6.82418(13) 


1.3 


-2.554998(13) 


-4.669821(60) 


-5.083112(75) 


-7.04020(14) 


1.4 


-2.643881(13) 


-4.827424(61) 


-5.253118(76) 


-7.27187(14) 


1.5 


-2.740369(13) 


-4.998013(62) 


-5.437176(78) 


-7.52257(14) 


1.6 


-2.846189(13) 


-5.184499(63) 


-5.638514(79) 


-7.79677(14) 


1.7 


-2.963693(14) 


-5.390845(65) 


-5.861540(81) 


-8.10050(16) 


1.8 


-3.096272(14) 


-5.622737(66) 


-6.112531(92) 


-8.44275(17) 


1.9 


-3.249289(14) 


-5.889063(76) 


-6.401622(95) 


-8.83741(17) 



Table 18. Numerical values of integrals I^^ and J^^ (Iwasaki). 
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-PL 




Ms 


unsmear 


MYPl 




0.1 


—2.592328(5) 


—5.035844(11) 


—4.947384(15) 


—6.309869(24) 


0.2 


—2.524721(4) 


—4.910064(9) 


—4.808736(12) 


—6.114771(21) 


0.3 


-2.460350(4) 


—4.790210(8) 


—4.677204(10) 


—5.930445(19) 


0.4 


—2.397926(4) 


—4.673907(7) 


—4.550068(9) 


—5.752916(17) 


0.5 


—2.336642(3) 


—4.559672(6) 


—4.425635(9) 


—5.579716(16) 


O.D 


—2.275903(3) 


—4.446413(6) 


—4.302662(8) 


CT A r\r\r\r"'7 / 1 a\ 

—5.409067(14) 


( \ r-f 

0.7 


—2.215218(3) 


—A.666226[b) 


—4.180126(7) 


— 5.239489(13) 


0.8 


— 2.154132(3) 


— 4.219272(5) 


—4.057093(7) 


—5.069649(13) 


0.9 


— 2.092214(3) 


—4.103759(5) 


—3.932667(7) 


—4.898273(12) 


1.0 


—2.029001(3) 


—3.985839(5) 


—3.805908(6) 


—4.724030(12) 


1.1 


— 1.963993(3) 


—3.864597(5) 


—3.675807(8) 


— 4.545491(12) 


1.2 


— 1.896617(3) 


— 3.738979(5) 


— 3.541177(8) 


— 4.360992(13) 


1.3 


— 1.826186(3) 


— 3.607730(5) 


— 3.400642(7) 


— 4.168585(15) 


1 A 
1.4 


— i. ( OiOOO^O ) 


— O.4:0VOOO[0 ) 


— O.ZoZoZO[v ) 


— o.yDoyuo(^io ) 


1.5 


-1.672630(3) 


-3.321845(5) 


-3.094636(7) 


-3.749859(15) 


1.6 


-1.587022(3) 


-3.162681(6) 


-2.924114(8) 


-3.516375(14) 


1.7 


-1.493120(3) 


-2.988297(6) 


-2.736979(8) 


-3.259815(15) 


1.8 


-1.388092(3) 


-2.793513(7) 


-2.527403(10) 


-2.971830(19) 


1.9 


-1.267425(4) 


-2.570075(8) 


-2.286043(13) 


-2.639040(25) 




Ms 


unsmear 


x' 

APE 


-PL 

HYPl 


HYP2 


0.1 


-1.615358(14) 


-3.028342(20) 


-3.226072(23) 


-4.411409(36) 


0.2 


-1.666709(15) 


-3.123267(21) 


-3.328015(24) 


-4.551286(37) 


0.3 


-1.717513(15) 


-3.217135(21) 


-3.428345(24) 


-4.688303(37) 


0.4 


-1.768504(15) 


-3.311288(22) 


-3.528584(24) 


-4.824649(38) 


0.5 


-1.820189(16) 


-3.406644(22) 


-3.629759(25) 


-4.961784(39) 


0.6 


-1.872985(16) 


-3.503956(22) 


-3.732706(25) 


-5.100881(39) 


0.7 


-1.927281(16) 


-3.603921(23) 


-3.838190(26) 


-5.243002(40) 


0.8 


-1.983471(16) 


-3.707237(23) 


-3.946969(26) 


-5.389194(40) 


0.9 


-2.041971(17) 


-3.814645(23) 


-4.059846(26) 


-5.540555(41) 


1.0 


-2.103248(17) 


-3.926964(24) 


-4.177706(27) 


-5.698292(42) 


1.1 


-2.167837(17) 


-4.045132(24) 


-4.301562(27) 


-5.863785(42) 


1.2 


-2.236373(17) 


-4.170254(24) 


-4.432607(28) 


-6.038655(43) 


1.3 


-2.309624(18) 


-4.303666(25) 


-4.572285(28) 


-6.224868(43) 


1.4 


-2.388547(18) 


-4.447019(25) 


-4.722389(29) 


-6.424869(44) 


1.5 


-2.474366(18) 


-4.602419(26) 


-4.885212(29) 


-6.641789(45) 


1.6 


-2.568694(19) 


-4.772629(26) 


-5.063771(30) 


-6.879766(45) 


1.7 


-2.673747(19) 


-4.961423(27) 


-5.262215(30) 


-7.144503(46) 


1.8 


-2.792740(20) 


-5.174249(27) 


-5.486561(31) 


-7.444311(47) 


1.9 


-2.930820(23) 


-5.419755(30) 


-5.746415(33) 


-7.792505(50) 



Table 19. Numerical values of integrals and J^} (Iwasaki). 
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-PL 




Ms 


unsmear 


hi YFl 


H YFz 


0.1 


10.096015(32) 


19.775670(48) 


18.955169(49) 


23.600213(75) 


0.2 


4.194416(21) 


8.271905(28) 


7.773551(32) 


9.484512(51) 


O.o 


2.176515(18) 


4.341016(27) 


o n /I onn /'o'7\ 

3.946890(27) 


4.646750(42) 


0.4 


1.125308(16) 


2.295156(25) 


1 nri 1 on/'on\ 

1.951130(30) 


2.118736(49) 


0.5 


0.456155(14) 


0.994411(20) 


0.678979(23) 


0.503499(35) 


O.D 


—0.027125(14) 


0.056320(19) 


—0.241092(22) 


—0.667694(33) 


0.7 


—0.410081(14) 


—0.685851(18) 


—0.971164(21) 


1 crnn /I /^n /'oo\ 

— 1.599469(32) 


0.8 


—0.737285(14) 


— 1.318921(19) 


— 1.595696(21) 


—2.398537(33) 


0.9 


— 1.035791(14) 


— 1.895458(20) 


—2.166002(23) 


—3.129847(35) 


1.0 


— 1.324912(15) 


—2.452964(21) 


—2.718689(24) 


—3.839838(37) 


1.1 


— 1.621273(16) 


—3.023576(23) 


— 3.285352(26) 


—4.568746(40) 


1.2 


— 1.942401(17) 


—3.641067(24) 


o onnon /i /o^\ 

—3.899304(27) 


—5.359136(42) 


1.3 


—2.310611(18) 


—4.348303(26) 


A /^non'^n /'on\ 

—4.602970(29) 


— 6.265350(45) 


1 4 


—9 75Q141 ('90'! 


— 5 9nQn95('98'l 


— 5 45Q55fi('31 


— 7 3fi8470('4Q'l 


1.5 


-3.344628(22) 


-6.331731(30) 


-6.576745(34) 


-8.806717(53) 


1.6 


-4.177540(24) 


-7.927926(33) 


-8.164523(38) 


-10.849676(59) 


1.7 


-5.512888(28) 


-10.485721(45) 


-10.707531(53) 


-14.119546(86) 


1.8 


-8.112468(35) 


-15.462945(55) 


-15.653132(65) 


-20.47423(11) 


1.9 


-15.779753(55) 


-30.137524(83) 


-30.226655(98) 


-39.18792(12) 



Table 20. Numerical values of integrals IC^ (Iwasaki). 
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Ms 


unsmear 




-PL 

MYPl 






0.1 


0.414157(5) 


0.491310(6) 


0.490776(6] 


0.537070(7] 




0.2 


0.426612(6] 


0.505961(6] 


0.505541(6] 


0.553354(7] 




0.3 


0.439797(6] 


r\ croi Acv (n\ 

0.521465(6] 


0.521169(6] 


0.570591(8] 




0.4 


0.453766(6] 


0.537884(7) 


0.537723(7] 


0.588849(8] 




0.5 


0.468581(6] 


0.555291(7) 


0.555275(7] 


0.608210(8] 




0.6 


0.484312(6] 


0.573767(7) 


0.573909(7] 


0.628763(9] 




0.7 


r\ v r\-\ r\ A A ( 'l\ 

0.501044(7] 


0.593408(7) 


0.593721(7] 


0.650617(9] 




0.8 


0.518870(7] 


0.614325(8) 


0.614823(8] 


0.673894(9] 




0.9 


0.537902(7] 


0.636644(8) 


0.637346(8] 


0.698737(10] 




1.0 


0.558268(7] 


0.660516(8) 


0.661439(8] 


0.725315(10] 




1.1 


0.580121(8] 


0.686114(8) 


0.687283(9] 


0.753824(10] 




1.2 


0.603641(8] 


0.713647(9] 


0.715086(9] 


0.784496(11] 




1.3 


0.629041(8] 


0.743359(9] 


0.745101(9] 


0.817612(11] 




1.4 


U.DODO /^y(^y ] 


[j. 1 1 oow[W ) 


u. / ( /Dzy^iu ] 


U.ooooUDi^iz ] 




1.5 


0.686569(9) 


0.810572(10) 


0.813038(10) 


0.892585(12) 




1.6 


0.719398(10) 


0.848874(11) 


0.851780(11) 


0.935354(13) 




1.7 


0.755548(10) 


0.891007(11) 


0.894421(11) 


0.982442(14) 




1.8 


0.795635(11) 


0.937670(12) 


0.941681(12) 


1.034650(15) 




1.9 


0.840457(12) 


0.989774(13) 


0.994493(13) 


1.093020(16) 




Ms 




unsmear 


APE 


-PL 

HYPl 


HYP2 


0.1 


-4 


956465(34) - 


-5.948112(53) 


-6.184557(61) 


-7.175375(98) 


0.2 


-5 


111478(35) - 


-6.133504(54) 


-6.377501(62) 


-7.39916(10) 


0.3 


-5 


260412(36) - 


-6.311567(55) 


-6.562570(63) 


-7.61340(10) 


0.4 


-5 


406144(36) - 


-6.485741(56) 


-6.743383(64) 


-7.82233(10) 


0.5 


-5 


550525(37) - 


-6.658231(57) 


-6.922257(65) 


-8.02870(10) 


0.6 


-5 


694973(38) - 


-6.830728(57) 


-7.100969(66) 


-8.23457(11) 


0.7 


-5 


840719(38) - 


-7.004694(58) 


-7.281045(67) 


-8.44173(11) 


0.8 


-5 


988924(39) - 


-7.181506(59) 


-7.463925(68) 


-8.65183(11) 


0.9 


-6 


140762(40) - 


-7.362553(60) 


-7.651053(69) 


-8.86656(11) 


1.0 


-6 


297486(40) - 


-7.549317(61) 


-7.843973(70) 


-9.08770(11) 


1.1 


-6 


460494(41) - 


-7.743445(62) 


-8.044400(71) 


-9.31722(11) 


1.2 


-6 


631408(41) - 


-7.946852(63) 


-8.254323(72) 


-9.55742(12) 


1.3 


-6 


812174(42) - 


-8.161828(64) 


-8.476125(73) 


-9.81103(12) 


1.4 


-7 


005202(43) - 


-8.391209(64) 


-8.712756(75) 


-10.08146(12) 


1.5 


-7 


213574(43) - 


-8.638618(66) 


-8.967989(76) 


-10.37305(12) 


1.6 


-7 


441377(44) - 


-8.908857(67) 


-9.246828(77) 


-10.69157(12) 


1.7 


-7 


694270(45) - 


-9.208571(68) 


-9.556204(78) 


-11.04503(12) 


1.8 


-7 


980554(46) - 


-9.547508(69) 


-9.906295(80) 


-11.44519(13) 


1.9 


-8 


313591(47) - 


-9.941352(71) 


-10.313452(89) 


-11.91096(14) 



Table 21. Numerical values of integrals I^^-^^ and J^^'^^ (DBW2). 
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A TDT7 


-PL 




Ms 


unsmear 


11 YPl 


riXrl 


0.1 


—5.079951(12) 


—6.202083(14) 


—6.210530(15) 


—6.914069(20) 


0.2 


—4.883963(11) 


— 5.966157(12) 


— 5.964620(13) 


—6.627483(17) 


0.3 


A '1f\ A A nf\( '\ cW 

—4.704460(10) 


— 5.749961(11) 


—5.739652(12) 


—6.365839(15) 


n A 

0.4 


—4.535904(9) 


— 5.546866(10) 


—5.528637(10) 


—6.120876(13) 


0.5 


—4.374955(8) 


—5.352871(9) 


—5.327356(9) 


—5.887615(12) 


O.D 


A A'1 ( 'l\ 

—4:.2l92i7[7) 


—5.165142(8) 


—5.132825(8) 


—5.662537(11) 


0.7 


—4.066914(6) 


— 4.y8i44z(7j 


A f\ A a f\r\ ( 'y\ 

—4.942690(7) 


— 5.442571(iUj 


0.8 


—3.916351(6) 


—4.799849(6) 


—4.754935(7) 


—5.226247(9) 


0.9 


—3.766078(5) 


—4.618585(6) 


—4.567695(6) 


—5.010482(8) 


1.0 


—3.614629(5) 


—4.435897(5) 


—4.379133(6) 


—4.793423(7) 


1.1 


o Anr\Ar'A/r'\ 

— 3.460464(5) 


—4.249928(5) 


—4.187317(5) 


—4.572813(7) 


1.2 


—3.301862(5) 


—4.058623(5) 


—3.990094(6) 


—4.346134(9) 




— 3.136830(5) 


—3.859588(5) 


—3.784955(6) 


A 1 1 r\ A '71 /'o \ 

— 4.1iU47i(oj 


1 A 
1.4 


— z.yDzoyu^o j 


— o.D4yo4y(^D j 


— o.ODooUo^D j 




1.5 


-2.776871(6) 


-3.425602(6) 


-3.337669(7) 


-3.596696(9) 


1.6 


-2.574495(6) 


-3.181723(7) 


-3.086179(7) 


-3.307664(9) 


1.7 


-2.349691(7) 


-2.910932(8) 


-2.806716(8) 


-2.986193(11) 


1.8 


-2.093217(9) 


-2.602158(10) 


-2.487666(10) 


-2.618670(13) 


1.9 


-1.789501(9) 


-2.236747(11) 


-2.109458(11) 


-2.182115(15) 




Ms 


unsmear 


x' 

APE 


-PL 

HYPl 


HYP2 


0.1 


-4.281982(49) 


-5.161542(50) 


-5.335295(51) 


-6.156654(56) 


0.2 


-4.413021(50) 


-5.319005(51) 


-5.498172(52) 


-6.344465(57) 


0.3 


-4.537621(51) 


-5.468743(52) 


-5.652733(53) 


-6.522188(58) 


0.4 


-4.658679(52) 


-5.614217(53) 


-5.802622(54) 


-6.694118(59) 


0.5 


-4.778024(52) 


-5.757616(54) 


-5.950136(55) 


-6.862948(60) 


0.6 


-4.897031(53) 


-5.900582(55) 


-6.096996(56) 


-7.030691(61) 


0.7 


-5.016865(54) 


-6.044503(56) 


-6.244650(56) 


-7.199031(62) 


0.8 


-5.138600(55) 


-6.190657(56) 


-6.394429(57) 


-7.369510(62) 


0.9 


-5.263302(56) 


-6.340312(57) 


-6.547649(58) 


-7.543645(63) 


1.0 


-5.392091(56) 


-6.494800(58) 


-6.705689(59) 


-7.723024(64) 


1.1 


-5.526209(57) 


-6.655589(59) 


-6.870071(60) 


-7.909392(65) 


1.2 


-5.667084(58) 


-6.824374(59) 


-7.042546(60) 


-8.104757(66) 


1.3 


-5.816428(58) 


-7.003175(60) 


-7.225207(61) 


-8.311515(67) 


1.4 


-5.976357(59) 


-7.194495(61) 


-7.420644(62) 


-8.532634(68) 


1.5 


-6.149592(60) 


-7.401542(62) 


-7.632181(63) 


-8.771926(68) 


1.6 


-6.339756(61) 


-7.628587(63) 


-7.864251(64) 


-9.034481(69) 


1.7 


-6.551909(62) 


-7.881590(64) 


-8.123045(65) 


-9.327421(71) 


1.8 


-6.793565(63) 


-8.169390(65) 


-8.417768(66) 


-9.661352(72) 


1.9 


-7.077040(65) 


-8.506461(67) 


-8.763516(68) 


-10.053716(74) 



Table 22. Numerical values of integrals J^^"^^ and J^}^'^^ (DBW2). 
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unsmear 


APE 


-FL 

HYPl 


HYP2 


0.1 


17.744304(92) 


21.765134(94) 


21.488573(95) 


23.51426(10) 


0.2 


6.622969(64) 


8.165288(66) 


7.942122(67) 


8.529364(73) 


0.3 


2.798492(54) 


3.489908(56) 


3.281417(57) 


3.368939(62) 


0.4 


0.792183(50) 


1.038274(52) 


0.834997(53) 


0.657026(56) 


0.5 


-0.494667(47) 


-0.533436(48) 


-0.735208(49) 


-1.085924(53) 


0.6 


-1.431073(46) 


-1.676488(47) 


-1.878567(48) 


-2.356828(52) 


0.7 


-2.178162(41) 


-2.587914(42) 


-2.791338(43) 


-3.372775(48) 


0.8 


-2.820196(46) 


-3.370709(48) 


-3.576167(48) 


-4.247368(53) 


0.9 


-3.408360(49) 


-4.087399(50) 


-4.295395(51) 


-5.049632(55) 


1.0 


-3.979486(51) 


-4.782965(53) 


-4.993953(53) 


-5.829417(58) 


1.1 


-4.565450(54) 


-5.496271(56) 


-5.710707(56) 


-6.629888(62) 


1.2 


-5.200017(57) 


-6.268436(59) 


-6.486849(60) 


-7.496864(65) 


1.3 


-5.926315(61) 


-7.151936(62) 


-7.375000(63) 


-8.488941(69) 


1.4 


-6.808610(65) 


-8.224926(67) 


-8.453584(68) 


-9.693491(74) 


1.5 


-7.956387(69) 


-9.620503(71) 


-9.856188(73) 


-11.259369(79) 


1.6 


-9.582978(76) 


-11.597969(78) 


-11.843078(79) 


-13.476581(87) 


1.7 


-12.180249(86) 


-14.755134(89) 


-15.014268(90) 


-17.013640(98) 


1.8 


-17.21586(10) 


-20.87572(11) 


-21.15993(11) 


-23.86488(12) 


1.9 


-32.01091(15) 


-38.85726(16) 


-39.20924(16) 


-43.97686(19) 



Table 23. Numerical values of integrals JC^^-^^ (DBW2). 
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^1 MP 


(^5 ) 




o 


unsmear 


APE 


HYPl 


HYP2 


0.1 


-8.01062(55) 


-3.52540(55) 


-3.73250(55) 


-1.21639(56) 


0.2 


-7.95272(55) 


-3.45640(55) 


-3.66320(55) 


-1.13986(55) 


0.3 


-7.90711(55) 


-3.39905(55) 


-3.60549(55) 


-1.07445(55) 


0.4 


-7.87113(55) 


-3.35062(55) 


-3.55666(55) 


-1.01740(55) 


0.5 


-7.84305(55) 


-3.30932(55) 


-3.51489(55) 


-0.96686(55) 


0.6 


-7.82239(55) 


-3.27461(55) 


-3.47965(55) 


-0.92224(55) 


0.7 


-7.80870(55) 


-3.24596(55) 


-3.45040(55) 


-0.88292(55) 


0.8 


-7.80192(55) 


-3.22322(55) 


-3.42696(55) 


-0.84868(55) 


0.9 


-7.80223(55) 


-3.20646(55) 


-3.40941(55) 


-0.81948(55) 


1.0 


-7.81014(55) 


-3.19608(55) 


-3.39811(55) 


-0.79562(55) 


1.1 


-7.82618(55) 


-3.19248(55) 


-3.39346(55) 


-0.77737(55) 


1.2 


-7.85132(55) 


-3.19645(55) 


-3.39623(55) 


-0.76535(55) 


1.3 


-7.88683(55) 


-3.20909(55) 


-3.40747(55) 


-0.76047(55) 


1.4 


-7.93440(55) 


-3.23186(55) 


-3.42862(55) 


-0.76396(56) 


1.5 


-7.99637(55) 


-3.26683(55) 


-3.46169(55) 


-0.77760(56) 


1.6 


-8.07600(55) 


-3.31694(55) 


-3.50957(55) 


-0.80398(56) 


1.7 


-8.17803(55) 


-3.38653(55) 


-3.57650(55) 


-0.84695(56) 


1.8 


-8.30970(55) 


-3.48232(55) 


-3.66912(56) 


-0.91268(56) 


1.9 


-8.48333(56) 


-3.61599(56) 


-3.79894(56) 


-1.01205(56) 








Zp 








unsmear 


APE 


HYPl 


HYP2 


0.1 


5.60787(65) 


-0.71311(65) 


-0.74216(65) 


-4.67527(66) 


0.2 


5.25969(57) 


-1.05003(57) 


-1.07660(57) 


-4.99868(58) 


0.3 


5.08451(38) 


-1.21647(39) 


-1.24146(39) 


-5.15550(41) 


0.4 


4.97796(33) 


-1.31639(34) 


-1.34047(34) 


-5.24900(36) 


0.5 


4.90762(30) 


-1.38199(30) 


-1.40588(31) 


-5.31108(33) 


0.6 


4.86000(27) 


-1.42668(27) 


-1.45098(28) 


-5.35496(31) 


0.7 


4.82793(30) 


-1.45754(31) 


-1.48291(31) 


-5.38771(34) 


0.8 


4.80736(32) 


-1.47880(32) 


-1.50572(33) 


-5.41339(35) 


0.9 


4.79552(27) 


-1.49313(28) 


-1.52225(28) 


-5.43489(31) 


1.0 


4.79009(29) 


-1.50297(29) 


-1.53494(30) 


-5.45473(33) 


1.1 


4.78979(30) 


-1.50975(31) 


-1.54521(31) 


-5.47447(34) 


1.2 


4.79266(39) 


-1.51560(40) 


-1.55525(40) 


-5.49648(42) 


1.3 


4.79754(41) 


-1.52188(42) 


-1.56650(42) 


-5.52243(45) 


1.4 


4.80128(44) 


-1.53202(45) 


-1.58246(46) 


-5.55611(48) 


1.5 


4.80292(40) 


-1.54732(40) 


-1.60455(41) 


-5.59937(44) 


1.6 


4.79352(43) 


-1.57711(44) 


-1.64232(44) 


-5.66225(47) 


1.7 


4.76193(47) 


-1.63330(48) 


-1.70774(48) 


-5.75744(51) 


1.8 


4.67607(65) 


-1.74875(66) 


-1.83405(66) 


-5.91919(68) 


1.9 


4.43275(88) 


-2.02799(88) 


-2.12627(88) 


-6.25418(90) 



Table 24. Numerical values of -SpV^^ and •''^^ (Plaquette). 



- 62 - 





unsmear 


V, [rna jNir 

APE 


HYPl 


HYP2 


0.1 


-36.2208(17) 


-24.6882(17) 


-26.4405(17) 


-22.1150(17) 


0.2 


-15.26247(82) 


-13.76112(82) 


-14.64156(82) 


-15.11317(82) 


0.3 


-8.07443(58) 


-9.98326(58) 


-10.57070(58) 


-12.67842(58) 


0.4 


-4.31899(47) 


-7.99166(47) 


-8.43052(47) 


-11.38912(47) 


0.5 


-1.92165(50) 


-6.70836(51) 


-7.05620(51) 


-10.55654(51) 


0.6 


-0.18470(39) 


-5.76989(39) 


-6.05514(39) 


-9.94787(39) 


0.7 


1.19577(33) 


-5.01790(33) 


-5.25640(33) 


-9.46192(33) 


0.8 


2.37853(32) 


-4.36917(32) 


-4.57049(32) 


-9.04562(33) 


0.9 


3.46160(34) 


-3.77143(34) 


-3.94128(34) 


-8.66515(34) 


1.0 


4.51340(40) 


-3.18863(40) 


-3.33040(40) 


-8.29815(40) 


1.1 


5.59437(39) 


-2.58791(39) 


-2.70310(39) 


-7.92398(40) 


1.2 


6.76817(48) 


-1.93436(48) 


-2.02281(48) 


-7.52110(48) 


1.3 


8.11661(56) 


-1.18244(56) 


-1.24216(56) 


-7.06146(56) 


1.4 


9.76147(51) 


-0.26402(51) 


-0.29050(51) 


-6.50347(52) 


1.5 


11.90951(73) 


0.93660(74) 


0.95197(74) 


-5.77686(74) 


1.6 


14.96610(67) 


2.64753(67) 


2.72114(67) 


-4.74221(68) 


1.7 


19.86546(91) 


5.39428(92) 


5.56044(92) 


-3.07893(92) 


1.8 


29.3973(10) 


10.7464(10) 


11.0929(10) 


0.1709(10) 


1.9 


57.5033(18) 


26.5554(18) 


27.4376(18) 


9.8075(18) 





unsmear 


APE 


HYPl 


HYP2 


0.1 


-26.1654(15) 


-11.2147(15) 


-11.9050(15) 


-3.5180(15) 


0.2 


-25.7609(15) 


-10.7732(15) 


-11.4625(15) 


-3.0514(15) 


0.3 


-25.4394(15) 


-10.4125(15) 


-11.1006(15) 


-2.6638(15) 


0.4 


-25.1696(15) 


-10.1013(15) 


-10.7881(15) 


-2.3239(15) 


0.5 


-24.9361(15) 


-9.8236(15) 


-10.5089(15) 


-2.0155(15) 


0.6 


-24.7317(15) 


-9.5725(15) 


-10.2559(15) 


-1.7312(15) 


0.7 


-24.5514(15) 


-9.3422(15) 


-10.0237(15) 


-1.4654(15) 


0.8 


-24.3918(15) 


-9.1294(15) 


-9.8085(15) 


-1.2143(15) 


0.9 


-24.2509(15) 


-8.9317(15) 


-9.6082(15) 


-0.9751(15) 


1.0 


-24.1279(15) 


-8.7477(15) 


-9.4212(15) 


-0.7462(15) 


1.1 


-24.0211(15) 


-8.5754(15) 


-9.2454(15) 


-0.5250(15) 


1.2 


-23.9319(15) 


-8.4157(15) 


-9.0816(15) 


-0.3120(15) 


1.3 


-23.8592(15) 


-8.2667(15) 


-8.9280(15) 


-0.1046(15) 


1.4 


-23.8035(15) 


-8.1283(15) 


-8.7842(15) 


0.0980(15) 


1.5 


-23.7656(15) 


-8.0004(15) 


-8.6500(15) 


0.2970(15) 


1.6 


-23.7460(15) 


-7.8824(15) 


-8.5245(15) 


0.4941(15) 


1.7 


-23.7446(15) 


-7.7729(15) 


-8.4061(15) 


0.6924(15) 


1.8 


-23.7590(15) 


-7.6678(15) 


-8.2904(15) 


0.8977(15) 


1.9 


-23.7759(15) 


-7.5515(15) 


-8.1613(15) 


1.1283(15) 



Table 25. Numerical values of ^^""^^^ and ^^Vm^ (Plaquette). 



- 63 - 





unsmear 




t )i/rMF\ 

rl Y Fi 


XJ\/TDO 

Ml rz 


0.1 


— 11.2157(13) 


1.4262(13) 


1.4843(13) 


9.3505(13) 


0.2 




— 10.5194(11) 


o 1 r\r\ 1 / 1 1 \ 

2.1001(11) 


2.1532(11) 


9.9974(12) 


0.3 




-10.16903(77) 


2.43293(78) 


2.48293(78) 


10.31099(82) 


0.4 




—9.95591(66) 


2.63277(67) 


o C) r\r\ A f n c>\ 

2.68094(68) 


10.49801(72) 


0.5 




—9.81524(59) 


2.76398(60) 


2.81175(61) 


10.62216(66) 


0.6 


—9.72000(54) 


2.85336(55) 


o f\r\'\ r\''7 / v r*\ 

2.90197(56) 


10.70992(62) 


0.7 


—9.65586(60) 


2.91508(62) 


2.96582(63) 


10.77542(68) 


0.5 




—9.61472(63) 


2.95760(64) 


3.01145(66) 


10.82678(71) 


0.9 


—9.59104(54) 


2.98626(56) 


3.04451(57) 


10.86978(63) 


1.0 


—9.58019(57) 


3.00594(59) 


3.06988(60) 


10.90945(66) 


1.1 


—9.57958(60) 


3.01951(61) 


3.09042(62) 


10.94894(68) 


1.2 




—9.58532(78) 


3.03119(79) 


3.11050(80) 


10.99296(85) 


1.3 




-9.59507(82) 


3.04376(83) 


3.13299(84) 


11.04485(89) 


1.4 




— y.DU/oD(oy ) 


o.UD4Uo(yU ) 


o.iD4yz(yi ) 


11 11 ooo/'nf^\ 
ii.iizzz(^yD ) 


1.5 




-9.60584(80) 


3.09464(81) 


3.20911(82) 


11.19874(87) 


1.6 


-9.58705(86) 


3.15423(88) 


3.28465(88) 


11.32450(94) 


1.7 


-9.52386(94) 


3.26659(96) 


3.41549(96) 


11.5149(10) 


1.5 




-9.3521(13) 


3.4975(13) 


3.6681(13) 


11.8384(14) 


1.9 


-8.8655(18) 


4.0560(18) 


4.2525(18) 


12.5084(18) 




^MF 




unsmear 


APE 


(MMF) 

HYPl 


HYP2 


0.1 




72.4416(34) 


49.3765(34) 


52.8811(34) 


44.2299(34) 


0.2 




30.5249(16) 


27.5222(16) 


29.2831(16) 


30.2263(16) 


0.3 




16.1489(12) 


19.9665(12) 


21.1414(12) 


25.3568(12) 


0.4 




8.63798(94) 


15.98332(94) 


16.86104(94) 


22.77823(95) 


0.5 




3.8433(10) 


13.4167(10) 


14.1124(10) 


21.1131(10) 


0.6 




0.36940(78) 


11.53978(78) 


12.11027(78) 


19.89574(78) 


0.7 




-2.39154(66) 


10.03579(66) 


10.51280(66) 


18.92384(67) 


0.8 




-4.75706(64) 


8.73834(64) 


9.14097(65) 


18.09124(66) 


0.9 




-6.92321(68) 


7.54286(68) 


7.88257(68) 


17.33029(69) 


1.0 




-9.02680(80) 


6.37726(80) 


6.66080(80) 


16.59630(81) 


1.1 




-11.18874(78) 


5.17582(78) 


5.40620(78) 


15.84796(79) 


1.2 




-13.53635(96) 


3.86872(96) 


4.04562(96) 


15.04219(97) 


1.3 




-16.2332(11) 


2.3649(11) 


2.4843(11) 


14.1229(11) 


1.4 




-19.5229(10) 


0.5280(10) 


0.5810(10) 


13.0069(10) 


1.5 




-23.8190(15) 


-1.8732(15) 


-1.9039(15) 


11.5537(15) 


1.6 




-29.9322(13) 


-5.2951(13) 


-5.4423(13) 


9.4844(14) 


1.7 




-39.7309(18) 


-10.7886(18) 


-11.1209(18) 


6.1579(18) 


1.8 




-58.7946(21) 


-21.4928(21) 


-22.1858(21) 


-0.3418(21) 


1.9 




-115.0066(36) 


-53.1107(36) 


-54.8752(36) 


-19.6151(36) 


Table 26. Numerical values of zj^^"^ 


VIF J -(ma)MF 

and 


(Plaquette). 



-64- 







^1 MP 


(A^5 ) 




o 


unsmear 


APE 


HYPl 


HYP2 


0.1 


-7.12281(55) 


-3.02368(55) 


-3.08112(55) 


-1.00130(56) 


0.2 


-7.06468(55) 


-2.95680(55) 


-3.01398(55) 


-0.92842(55) 


0.3 


-7.01857(55) 


-2.90144(55) 


-2.95831(55) 


-0.86663(55) 


0.4 


-6.98178(55) 


-2.85484(55) 


-2.91137(55) 


-0.81317(55) 


0.5 


-6.95256(55) 


-2.81523(55) 


-2.87136(55) 


-0.76620(55) 


0.6 


-6.93034(55) 


-2.78195(55) 


-2.83763(55) 


-0.72504(55) 


0.7 


-6.91468(55) 


-2.75452(55) 


-2.80969(55) 


-0.68914(55) 


0.8 


-6.90544(55) 


-2.73273(55) 


-2.78733(55) 


-0.65824(55) 


0.9 


-6.90272(55) 


-2.71660(55) 


-2.77055(55) 


-0.63228(55) 


1.0 


-6.90694(55) 


-2.70647(55) 


-2.75968(55) 


-0.61150(55) 


1.1 


-6.91866(55) 


-2.70278(55) 


-2.75514(55) 


-0.59626(55) 


1.2 


-6.93857(55) 


-2.70611(55) 


-2.75750(55) 


-0.58700(55) 


1.3 


-6.96786(55) 


-2.71749(55) 


-2.76778(55) 


-0.58463(55) 


1.4 


-7.00808(55) 


-2.73832(55) 


-2.78732(55) 


-0.59035(55) 


1.5 


-7.06140(55) 


-2.77054(55) 


-2.81806(55) 


-0.60589(55) 


1.6 


-7.13077(55) 


-2.81687(55) 


-2.86264(55) 


-0.63369(56) 


1.7 


-7.22060(55) 


-2.88140(55) 


-2.92511(55) 


-0.67750(56) 


1.8 


-7.33770(55) 


-2.97058(55) 


-3.01182(55) 


-0.74330(56) 


1.9 


-7.49384(56) 


-3.09565(56) 


-3.13392(56) 


-0.84177(56) 








Zp 








unsmear 


APE 


HYPl 


HYP2 


0.1 


4.10388(65) 


-1.29859(65) 


-1.40715(65) 


-4.90191(66) 


0.2 


3.75809(57) 


-1.63258(57) 


-1.73850(57) 


-5.22164(58) 


0.3 


3.58447(38) 


-1.79682(39) 


-1.90095(39) 


-5.37544(40) 


0.4 


3.47878(33) 


-1.89512(34) 


-1.99814(34) 


-5.46624(36) 


0.5 


3.40868(30) 


-1.95959(30) 


-2.06212(30) 


-5.52588(32) 


0.6 


3.36069(27) 


-2.00360(27) 


-2.10618(28) 


-5.56753(30) 


0.7 


3.32764(30) 


-2.03424(31) 


-2.13742(31) 


-5.59817(33) 


0.8 


3.30546(32) 


-2.05557(32) 


-2.15986(32) 


-5.62181(34) 


0.9 


3.29137(27) 


-2.07041(28) 


-2.17632(28) 


-5.64126(30) 


1.0 


3.28300(29) 


-2.08114(29) 


-2.18921(30) 


-5.65888(33) 


1.1 


3.27901(30) 


-2.08918(30) 


-2.19994(31) 


-5.67634(34) 


1.2 


3.27738(39) 


-2.09666(40) 


-2.21068(40) 


-5.69583(42) 


1.3 


3.27684(41) 


-2.10497(42) 


-2.22280(42) 


-5.71890(44) 


1.4 


3.27415(44) 


-2.11752(45) 


-2.23988(45) 


-5.74920(48) 


1.5 


3.26818(40) 


-2.13563(40) 


-2.26326(41) 


-5.78837(43) 


1.6 


3.24980(43) 


-2.16872(44) 


-2.30243(44) 


-5.84623(46) 


1.7 


3.20756(47) 


-2.22855(48) 


-2.36927(48) 


-5.93504(50) 


1.8 


3.10901(65) 


-2.34795(66) 


-2.49690(66) 


-6.08847(68) 


1.9 


2.85031(88) 


-2.63160(88) 


-2.79006(88) 


-6.41215(90) 



Table 27. Numerical values of ^pVm^^ and ^^^"^^^ (Symanzik). 



- 65 - 





unsmear 


V, [rna jNir 

APE 


HYPl 


HYP2 


0.1 


-33.8405(17) 


-24.5697(17) 


-25.9817(17) 


-22.5181(17) 


0.2 


-15.07855(82) 


-14.09644(82) 


-14.86596(82) 


-15.48564(82) 


0.3 


-8.63821(58) 


-10.47437(58) 


-11.02831(58) 


-13.04088(58) 


0.4 


-5.27038(47) 


-8.56457(47) 


-9.00952(47) 


-11.74672(47) 


0.5 


-3.11868(50) 


-7.33420(51) 


-7.71267(51) 


-10.91136(51) 


0.6 


-1.55854(39) 


-6.43476(39) 


-6.76778(39) 


-10.30093(39) 


0.7 


-0.31796(33) 


-5.71452(33) 


-6.01390(33) 


-9.81378(33) 


0.8 


0.74523(32) 


-5.09384(32) 


-5.36668(32) 


-9.39656(33) 


0.9 


1.71898(34) 


-4.52252(34) 


-4.77320(34) 


-9.01540(34) 


1.0 


2.66450(40) 


-3.96616(40) 


-4.19732(40) 


-8.64775(40) 


1.1 


3.63606(39) 


-3.39337(39) 


-3.60632(39) 


-8.27290(40) 


1.2 


4.69082(48) 


-2.77085(48) 


-2.96576(48) 


-7.86925(48) 


1.3 


5.90228(56) 


-2.05522(56) 


-2.23099(56) 


-7.40867(56) 


1.4 


7.37992(51) 


-1.18165(51) 


-1.33552(51) 


-6.84943(51) 


1.5 


9.30959(73) 


-0.04012(74) 


-0.16666(74) 


-6.12121(74) 


1.6 


12.05590(67) 


1.58642(67) 


1.49775(67) 


-5.08440(68) 


1.7 


16.45906(91) 


4.19784(92) 


4.16922(92) 


-3.41805(92) 


1.8 


25.0283(10) 


9.2871(10) 


9.3755(10) 


-0.1634(10) 


1.9 


50.3053(18) 


24.3234(18) 


24.7599(18) 


9.4833(18) 



^MF 


unsmear 


.(1)MF 
APE 


HYPl 


HYP2 


0.1 


-23.6446(15) 


-9.9808(15) 


-10.1723(15) 


-3.2396(15) 


0.2 


-23.2451(15) 


-9.5522(15) 


-9.7428(15) 


-2.7909(15) 


0.3 


-22.9281(15) 


-9.2043(15) 


-9.3939(15) 


-2.4216(15) 


0.4 


-22.6623(15) 


-8.9059(15) 


-9.0943(15) 


-2.1003(15) 


0.5 


-22.4322(15) 


-8.6411(15) 


-8.8282(15) 


-1.8110(15) 


0.6 


-22.2305(15) 


-8.4025(15) 


-8.5881(15) 


-1.5461(15) 


0.7 


-22.0520(15) 


-8.1848(15) 


-8.3688(15) 


-1.3003(15) 


0.8 


-21.8935(15) 


-7.9845(15) 


-8.1664(15) 


-1.0695(15) 


0.9 


-21.7527(15) 


-7.7990(15) 


-7.9788(15) 


-0.8513(15) 


1.0 


-21.6286(15) 


-7.6270(15) 


-7.8044(15) 


-0.6438(15) 


1.1 


-21.5196(15) 


-7.4667(15) 


-7.6413(15) 


-0.4450(15) 


1.2 


-21.4268(15) 


-7.3185(15) 


-7.4899(15) 


-0.2549(15) 


1.3 


-21.3485(15) 


-7.1806(15) 


-7.3482(15) 


-0.0711(15) 


1.4 


-21.2855(15) 


-7.0530(15) 


-7.2163(15) 


0.1069(15) 


1.5 


-21.2380(15) 


-6.9351(15) 


-7.0935(15) 


0.2804(15) 


1.6 


-21.2061(15) 


-6.8264(15) 


-6.9789(15) 


0.4509(15) 


1.7 


-21.1893(15) 


-6.7253(15) 


-6.8710(15) 


0.6211(15) 


1.8 


-21.1843(15) 


-6.6273(15) 


-6.7647(15) 


0.7970(15) 


1.9 


-21.1773(15) 


-6.5167(15) 


-6.6443(15) 


0.9963(15) 



Table 28. Numerical values of ^j™")'^^ and z^V^^^ (Symanzik). 
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unsmear 




I )i /fMF\ 

XI Y Jr^i 


riXr I 


0. 


1 




—8.2078(13) 


2.5972(13) 


2.8143(13) 


9.8038(13) 


0. 


2 




— 7.5162(11) 


3.2652(11) 


3.4770(11) 


10.4433(12) 


0. 


3 




— 7.16894(77) 


3.59365(78) 


3.80191(78) 


10.75087(81) 


0. 


4 




—6.95757(66) 


3.79024(67) 


3.99628(68) 


10.93248(71) 


0. 


5 




—6.81737(59) 


3.91918(60) 


4.12423(61) 


11.05176(64) 


0. 


5 




—6.72138(54) 


4.00719(55) 


4.21237(55) 


11.13506(60) 


0. 


7 




—6.65528(60) 


4.06848(62) 


4.27484(62) 


11.19634(66) 


0. 


3 




—6.61093(63) 


4.11114(64) 


4.31972(65) 


11.24362(69) 


0. 


9 




—6.58274(54) 


4.14082(55) 


4.35264(56) 


11.28252(61) 


1. 


D 




r* c r* f\r\ ( ''7\ 

—6.56600(57) 


4.16228(58) 


4.37842(59) 


1 1 01 'I'lO I r'v\ 

11.31776(65) 


1. 


1 




—6.55802(60) 


4.17837(61) 


4.39989(62) 


11.35269(68) 


1. 


2 




—6.55475(78) 


4.19333(79) 


4.42136(80) 


1 1 ^C\'\ 1 o a\ 

11.39167(84) 


1. 


3 




—6.55368(82) 


4.20994(83) 


4.44oDU(o4j 


11.43779(89) 


1. 


4 




— D.o4ooU(^oy j 


A OQ^^/l/'^^^ 
4.zooU4(yu j 




ii.4yo4U(yu j 


1. 


5 




-6.53636(80) 


4.27127(81) 


4.52652(82) 


11.57674(87) 


1. 


S 




-6.49959(86) 


4.33744(87) 


4.60486(88) 


11.69246(93) 


1. 


7 




-6.41512(94) 


4.45709(95) 


4.73855(96) 


11.8701(10) 


1. 


B 




-6.2180(13) 


4.6959(13) 


4.9938(13) 


12.1769(14) 


1. 


9 




-5.7006(18) 


5.2632(18) 


5.5801(18) 


12.8243(18) 




^MF 






unsmear 


APE 


(MMF) 

HYPl 


HYP2 


0.1 




67.6810(34) 


49.1395(34) 


51.9635(34) 


45.0361(34) 


0.2 




30.1571(16) 


28.1929(16) 


29.7319(16) 


30.9713(16) 


0.3 




17.2764(12) 


20.9487(12) 


22.0566(12) 


26.0818(12) 


0.4 




10.54075(94) 


17.12915(94) 


18.01905(94) 


23.49345(95) 


0.5 






6.2373(10) 


14.6684(10) 


15.4253(10) 


21.8227(10) 


0.6 




3.11708(78) 


12.86952(78) 


13.53557(78) 


20.60186(78) 


0.7 




0.63592(66) 


11.42903(66) 


12.02780(66) 


19.62755(67) 


0.8 






-1.49046(64) 


10.18767(64) 


10.73336(65) 


18.79312(65) 


0.9 






-3.43796(68) 


9.04503(68) 


9.54640(68) 


18.03079(69) 


1.0 






-5.32900(80) 


7.93232(80) 


8.39465(80) 


17.29550(80) 


1.1 






-7.27213(78) 


6.78673(78) 


7.21265(78) 


16.54580(79) 


1.2 






-9.38163(96) 


5.54170(96) 


5.93152(96) 


15.73850(97) 


1.3 






-11.8046(11) 


4.1104(11) 


4.4620(11) 


14.8173(11) 


1.4 






-14.7598(10) 


2.3633(10) 


2.6710(10) 


13.6989(10) 


1.5 






-18.6192(15) 


0.0802(15) 


0.3333(15) 


12.2424(15) 


1.6 






-24.1118(13) 


-3.1728(13) 


-2.9955(13) 


10.1688(14) 


1.7 






-32.9181(18) 


-8.3957(18) 


-8.3384(18) 


6.8361(18) 


1.8 






-50.0566(21) 


-18.5742(21) 


-18.7510(21) 


0.3268(21) 


1.9 






100.6107(36) 


-48.6468(36) 


-49.5198(36) 


-18.9667(36) 



Table 29. Numerical values of and z\^°'' (Symanzik). 
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o 


unsmear 


^1 MP 
APE 


(A^5 ) 

HYPl 


HYP2 


0.1 


-5.05539(55) 


-1.86326(55) 


-1.78838(55) 


-0.27208(55) 


0.2 


-4.99901(55) 


-1.80139(55) 


-1.72631(55) 


-0.20633(55) 


0.3 


-4.95413(55) 


-1.75071(55) 


-1.67539(55) 


-0.15151(55) 


0.4 


-4.91804(55) 


-1.70847(55) 


-1.63289(55) 


-0.10486(55) 


0.5 


-4.88882(55) 


-1.67273(55) 


-1.59687(55) 


-0.06441(55) 


0.6 


-4.86595(55) 


-1.64294(55) 


-1.56676(55) 


-0.02959(55) 


0.7 


-4.84883(55) 


-1.61847(55) 


-1.54192(55) 


0.00029(55) 


0.8 


-4.83727(55) 


-1.59908(55) 


-1.52214(55) 


0.02547(55) 


0.9 


-4.83124(55) 


-1.58469(55) 


-1.50730(55) 


0.04609(55) 


1.0 


-4.83100(55) 


-1.57552(55) 


-1.49764(55) 


0.06197(55) 


1.1 


-4.83692(55) 


-1.57186(55) 


-1.49342(55) 


0.07289(55) 


1.2 


-4.84962(55) 


-1.57427(55) 


-1.49520(55) 


0.07836(55) 


1.3 


-4.86995(55) 


-1.58352(55) 


-1.50373(55) 


0.07771(55) 


1.4 


-4.89923(55) 


-1.60081(55) 


-1.52020(55) 


0.06982(55) 


1.5 


-4.93920(55) 


-1.62777(55) 


-1.54622(55) 


0.05320(55) 


1.6 


-4.99248(55) 


-1.66686(55) 


-1.58423(55) 


0.02554(55) 


1.7 


-5.06279(55) 


-1.72164(55) 


-1.63773(55) 


-0.01648(55) 


1.8 


-5.15624(55) 


-1.79796(55) 


-1.71255(55) 


-0.07847(56) 


1.9 


-5.28359(56) 


-1.90631(56) 


-1.81908(56) 


-0.17056(56) 






unsmear 


APE 


HYPl 


HYP2 


0.1 


1.24050(65) 


-2.69648(65) 


-2.90181(65) 


-5.65092(66) 


0.2 


0.90312(57) 


-3.02203(57) 


-3.22460(57) 


-5.96189(58) 


0.3 


0.73579(38) 


-3.17980(39) 


-3.38039(39) 


-6.10856(40) 


0.4 


0.63469(33) 


-3.27313(34) 


-3.47235(34) 


-6.19351(35) 


0.5 


0.56772(30) 


-3.33388(30) 


-3.53226(30) 


-6.24818(32) 


0.6 


0.52153(27) 


-3.37522(27) 


-3.57321(28) 


-6.28551(29) 


0.7 


0.48903(30) 


-3.40414(31) 


-3.60217(31) 


-6.31231(33) 


0.8 


0.46622(32) 


-3.42459(32) 


-3.62305(32) 


-6.33241(34) 


0.9 


0.45030(27) 


-3.43931(28) 


-3.63858(28) 


-6.34849(30) 


1.0 


0.43892(29) 


-3.45064(29) 


-3.65109(30) 


-6.36284(31) 


1.1 


0.43073(30) 


-3.45994(30) 


-3.66195(31) 


-6.37682(33) 


1.2 


0.42364(39) 


-3.46929(40) 


-3.67321(40) 


-6.39251(41) 


1.3 


0.41635(41) 


-3.48001(42) 


-3.68622(42) 


-6.41125(43) 


1.4 


0.40551(44) 


-3.49548(45) 


-3.70437(45) 


-6.43649(46) 


1.5 


0.38992(40) 


-3.51694(40) 


-3.72889(41) 


-6.46951(42) 


1.6 


0.36031(43) 


-3.55366(44) 


-3.76910(44) 


-6.51963(45) 


1.7 


0.30512(47) 


-3.61721(48) 


-3.83659(48) 


-6.59838(50) 


1.8 


0.19171(65) 


-3.74018(66) 


-3.96386(66) 


-6.73851(68) 


1.9 


-0.08397(88) 


-4.02639(88) 


-4.25492(88) 


-7.04370(90) 



Table 30. Numerical values of ^pVm^^ and (Iwasaki). 
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unsmear 


V, [rna jNir 

APE 


HYPl 


HYP2 


0.1 


-30.3325(17) 


-24.5093(17) 


-25.4391(17) 


-23.3419(17) 


0.2 


-15.25949(82) 


-15.02390(82) 


-15.62567(82) 


-16.44402(82) 


0.3 


-10.07578(58) 


-11.74076(58) 


-12.23309(58) 


-14.04642(58) 


0.4 


-7.36011(47) 


-10.00903(47) 


-10.44651(47) 


-12.77782(47) 


0.5 


-5.62233(51) 


-8.89356(51) 


-9.29807(51) 


-11.95965(51) 


0.6 


-4.36071(39) 


-8.07875(39) 


-8.46116(39) 


-11.36234(39) 


0.7 


-3.35681(33) 


-7.42722(33) 


-7.79371(33) 


-10.88619(33) 


0.8 


-2.49636(32) 


-6.86690(32) 


-7.22123(32) 


-10.47885(33) 


0.9 


-1.70837(34) 


-6.35226(34) 


-6.69691(34) 


-10.10707(34) 


1.0 


-0.94376(40) 


-5.85237(40) 


-6.18890(40) 


-9.74876(40) 


1.1 


-0.15873(39) 


-5.33894(39) 


-5.66835(39) 


-9.38365(40) 


1.2 


0.69278(48) 


-4.78213(48) 


-5.10492(48) 


-8.99061(48) 


1.3 


1.67019(56) 


-4.14308(56) 


-4.45928(56) 


-8.54219(56) 


1.4 


2.86190(51) 


-3.36392(51) 


-3.67300(51) 


-7.99778(51) 


1.5 


4.41788(73) 


-2.34649(74) 


-2.64709(74) 


-7.28886(74) 


1.6 


6.63285(67) 


-0.89713(67) 


-1.18631(67) 


-6.27972(68) 


1.7 


10.18570(92) 


1.43002(92) 


1.15873(92) 


-4.65841(92) 


1.8 


17.1040(10) 


5.9666(10) 


5.7302(10) 


-1.4931(10) 


1.9 


37.5253(18) 


19.3760(18) 


19.2442(18) 


7.8839(18) 



^MF 


unsmear 


.(1)MF 
APE 


HYPl 


HYP2 


0.1 


-17.5944(15) 


-6.9540(15) 


-6.7044(15) 


-1.6501(15) 


0.2 


-17.2129(15) 


-6.5541(15) 


-6.3038(15) 


-1.2373(15) 


0.3 


-16.9127(15) 


-6.2346(15) 


-5.9836(15) 


-0.9040(15) 


0.4 


-16.6629(15) 


-5.9643(15) 


-5.7124(15) 


-0.6189(15) 


0.5 


-16.4472(15) 


-5.7269(15) 


-5.4740(15) 


-0.3659(15) 


0.6 


-16.2589(15) 


-5.5155(15) 


-5.2616(15) 


-0.1377(15) 


0.7 


-16.0924(15) 


-5.3245(15) 


-5.0693(15) 


0.0714(15) 


0.8 


-15.9442(15) 


-5.1502(15) 


-4.8937(15) 


0.2650(15) 


0.9 


-15.8120(15) 


-4.9902(15) 


-4.7322(15) 


0.4457(15) 


1.0 


-15.6945(15) 


-4.8429(15) 


-4.5833(15) 


0.6154(15) 


1.1 


-15.5898(15) 


-4.7063(15) 


-4.4448(15) 


0.7762(15) 


1.2 


-15.4988(15) 


-4.5810(15) 


-4.3174(15) 


0.9278(15) 


1.3 


-15.4196(15) 


-4.4648(15) 


-4.1988(15) 


1.0726(15) 


1.4 


-15.3522(15) 


-4.3575(15) 


-4.0888(15) 


1.2113(15) 


1.5 


-15.2967(15) 


-4.2585(15) 


-3.9867(15) 


1.3447(15) 


1.6 


-15.2525(15) 


-4.1671(15) 


-3.8916(15) 


1.4743(15) 


1.7 


-15.2180(15) 


-4.0808(15) 


-3.8012(15) 


1.6030(15) 


1.8 


-15.1894(15) 


-3.9952(15) 


-3.7104(15) 


1.7365(15) 


1.9 


-15.1513(15) 


-3.8937(15) 


-3.6030(15) 


1.8921(15) 


Table 31. Numerical values of ^j™")'^^ aj^^^ ^L(i?)^ 


(Iwasaki). 
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, (pa) MF 



unsmear 



APE 



(Mr) 



HYPl 



HYP2 



0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 
1.0 
1.1 
1.2 
1.3 
1.4 
1.5 
1.6 
1.7 
1.8 
1.9 



-2.4810(13) 
-1.8062(11) 
-1.47157(77) 
-1.26939(66) 
-1.13544(59) 
-1.04305(54) 
-0.97807(60) 
-0.93244(63) 
-0.90060(54) 
-0.87785(57) 
-0.86146(60) 
-0.84728(78) 
-0.83270(82) 
-0.81103(89) 
-0.77984(80) 
-0.72062(86) 
-0.61024(94) 
-0.3834(13) 
0.1679(18) 



5.3930(13) 
6.0441(11) 
6.35960(78) 
6.54626(67) 
6.66776(60) 
6.75044(55) 
6.80828(61) 
6.84918(64) 
6.87862(55) 
6.90128(58) 
6.91989(61) 
6.93859(79) 
6.96003(83) 
6.99097(90) 
7.03388(81) 
7.10731(87) 
7.23442(95) 
7.4804(13) 
8.0528(18) 



5.8036(13) 
6.4492(11) 
6.76078(78) 
6.94471(68) 
7.06451(60) 
7.14643(55) 
7.20434(62) 
7.24610(65) 
7.27717(56) 
7.30219(59) 
7.32390(62) 
7.34642(79) 
7.37243(84) 
7.40873(90) 
7.45778(81) 
7.53821(88) 
7.67317(96) 
7.9277(13) 
8.5098(18) 



11.3018(13) 
11.9238(12) 
12.21711(80) 
12.38702(70) 
12.49636(64) 
12.57103(59) 
12.62462(65) 
12.66482(68) 
12.69697(60) 
12.72568(63) 
12.75365(65) 
12.78502(82) 
12.82250(87) 
12.87299(93) 
12.93901(85) 
13.03925(91) 
13.19676(100) 
13.4770(14) 
14.0874(18) 



^MF 






(Mr) 






unsmear 


APE 


HYPl 


HYP2 


0.1 


60.6649(34) 


49.0186(34) 


50.8782(34) 


46.6837(34) 


0.2 


30.5190(16) 


30.0478(16) 


31.2513(16) 


32.8880(16) 


0.3 


20.1516(12) 


23.4815(12) 


24.4662(12) 


28.0929(12) 


0.4 


14.72022(94) 


20.01805(94) 


20.89302(94) 


25.55564(95) 


0.5 


11.2447(10) 


17.7871(10) 


18.5961(10) 


23.9193(10) 


0.6 


8.72142(78) 


16.15749(78) 


16.92232(78) 


22.72468(78) 


0.7 


6.71362(66) 


14.85445(66) 


15.58742(66) 


21.77238(67) 


0.8 


4.99272(64) 


13.73380(65) 


14.44246(65) 


20.95770(65) 


0.9 


3.41675(68) 


12.70452(68) 


13.39382(68) 


20.21414(69) 


1.0 


1.88752(80) 


11.70473(80) 


12.37780(80) 


19.49752(80) 


1.1 


0.31747(78) 


10.67787(78) 


11.33670(78) 


18.76730(79) 


1.2 


-1.38556(96) 


9.56426(96) 


10.20983(96) 


17.98122(97) 


1.3 


-3.3404(11) 


8.2862(11) 


8.9186(11) 


17.0844(11) 


1.4 


-5.7238(10) 


6.7278(10) 


7.3460(10) 


15.9956(10) 


1.5 


-8.8358(15) 


4.6930(15) 


5.2942(15) 


14.5777(15) 


1.6 


-13.2657(13) 


1.7943(13) 


2.3726(13) 


12.5594(14) 


1.7 


-20.3714(18) 


-2.8600(18) 


-2.3175(18) 


9.3168(18) 


1.8 


-34.2080(21) 


-11.9332(21) 


-11.4604(21) 


2.9861(21) 


1.9 


-75.0506(36) 


-38.7519(36) 


-38.4884(36) 


-15.7677(36) 



Table 32. Numerical values of and z ™ (Iwasaki). 
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,(1)MP 



unsmear 



APE 



HYPl 



HYP2 



0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 
1.0 
1.1 
1.2 
1.3 
1.4 
1.5 
1.6 
1.7 
1.8 
1.9 



-1.42622(55) 
-1.37916(55) 
-1.34243(55) 
-1.31323(55) 
-1.28963(55) 
-1.27099(55) 
-1.25662(55) 
-1.24612(55) 
-1.23944(55) 
-1.23663(55) 
-1.23782(55) 
-1.24324(55) 
-1.25355(55) 
-1.26953(55) 
-1.29231(55) 
-1.32377(55) 
-1.36665(55) 
-1.42558(55) 
-1.50938(56) 



0.28399(55) 
0.33324(55) 
0.37230(55) 
0.40394(55) 
0.43014(55) 
0.45152(55) 
0.46880(55) 
0.48239(55) 
0.49236(55) 
0.49867(55) 
0.50122(55) 
0.49982(55) 
0.49382(55) 
0.48249(55) 
0.46474(55) 
0.43876(55) 
0.40186(55) 
0.34951(55) 
0.27299(56) 



0.40401(55) 
0.45338(55) 
0.49256(55) 
0.52434(55) 
0.55068(55) 
0.57222(55) 
0.58967(55) 
0.60345(55) 
0.61361(55) 
0.62015(55) 
0.62295(55) 
0.62182(55) 
0.61612(55) 
0.60513(55) 
0.58777(55) 
0.56222(55) 
0.52584(55) 
0.47408(55) 
0.39826(56) 



1.27331(55) 
1.32419(55) 
1.36498(55) 
1.39847(55) 
1.42662(55) 
1.45008(55) 
1.46957(55) 
1.48552(55) 
1.49801(55) 
1.50703(55) 
1.51249(55) 
1.51423(55) 
1.51164(55) 
1.50401(55) 
1.49032(55) 
1.46880(55) 
1.43686(55) 
1.39005(55) 
1.31980(56) 





unsmear 


APE HYPl 


HYP2 


0.1 


-4.42855(65) 


-6.54233(65) 


-6.78722(65) 


-8.48158(65) 


0.2 


-4.73575(57) 


-6.83997(57) 


-7.08243(57) 


-8.76695(58) 


0.3 


-4.88036(39) 


-6.97702(39) 


-7.21771(39) 


-8.89472(40) 


0.4 


-4.96409(33) 


-7.05465(34) 


-7.29406(34) 


-8.96525(35) 


0.5 


-5.01786(30) 


-7.10348(30) 


-7.34199(30) 


-9.00870(31) 


0.6 


-5.05422(27) 


-7.13587(27) 


-7.37380(28) 


-9.03711(29) 


0.7 


-5.07970(30) 


-7.15820(31) 


-7.39580(31) 


-9.05667(32) 


0.8 


-5.09792(32) 


-7.17399(32) 


-7.41150(32) 


-9.07072(34) 


0.9 


-5.11132(27) 


-7.18562(28) 


-7.42323(28) 


-9.08152(29) 


1.0 


-5.12202(29) 


-7.19512(29) 


-7.43301(29) 


-9.09103(31) 


1.1 


-5.13117(30) 


-7.20358(30) 


-7.44192(31) 


-9.10024(32) 


1.2 


-5.14061(39) 


-7.21282(40) 


-7.45176(40) 


-9.11089(41) 


1.3 


-5.15147(41) 


-7.22388(42) 


-7.46354(42) 


-9.12396(43) 


1.4 


-5.16682(45) 


-7.23978(45) 


-7.48029(45) 


-9.14240(46) 


1.5 


-5.18752(40) 


-7.26130(40) 


-7.50274(41) 


-9.16682(42) 


1.6 


-5.22235(43) 


-7.29706(44) 


-7.53948(44) 


-9.20571(45) 


1.7 


-5.28203(47) 


-7.35757(48) 


-7.60099(48) 


-9.26929(49) 


1.8 


-5.39770(66) 


-7.47360(66) 


-7.71789(66) 


-9.38779(67) 


1.9 


-5.67034(88) 


-7.74535(88) 


-7.99016(88) 


-9.66033(89) 




Table 33. Numerical values of z 


(1)MF , 4pa)MF 


(DBW2). 
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unsmear 


.(rna)MF^^MF) 

APE HYPl 


HYP2 


0.1 


-27.8384(17) 


-25.8193(17) 


-26.2781(17) 


-25.7773(17) 


0.2 


-17.93649(82) - 


-18.38235(82) 


-18.78314(82) 


-19.70506(82) 


0.3 


-14.51802(58) - 


-15.80322(58) 


-16.18540(58) 


-17.59352(59) 


0.4 


-12.72139(48) - 


-14.44180(48) 


-14.81525(48) 


-16.47696(48) 


0.5 


-11.56930(51) - 


-13.56558(51) 


-13.93435(51) 


-15.75814(51) 


0.6 


-10.73205(39) - 


-12.92691(39) 


-13.29309(39) 


-15.23475(40) 


0.7 


-10.06617(34) - 


-12.41809(34) 


-12.78291(34) 


-14.81891(34) 


0.8 


-9.49662(33) - 


-11.98269(33) 


-12.34700(33) 


-14.46460(33) 


0.9 


-8 


97614(35) - 


-11.58469(35) 


-11.94914(35) 


-14.14216(35) 


1.0 


-8 


47281(40) - 


-11.20026(41) 


-11.56538(41) 


-13.83246(41) 


1.1 


-7.95775(40) - 


-10.80742(40) 


-11.17373(40) 


-13.51772(40) 


1.2 


-7.40079(49) - 


-10.38326(49) 


-10.75132(49) 


-13.17958(49) 


1.3 


-6.76296(57) 


-9.89809(57) 


-10.26855(57) 


-12.79431(57) 


1.4 


-5.98639(52) 


-9.30780(52) 


-9.68157(52) 


-12.32687(52) 


1.5 


-4.97334(74) 


-8.53814(74) 


-8.91653(74) 


-11.71848(74) 


1.6 


-3.53113(68) 


-7.44218(68) 


-7.82741(68) 


-10.85262(68) 


1.7 


-1.21639(92) 


-5.68220(92) 


-6.07857(92) 


-9.46180(92) 


1.8 


3.2947(10) 


-2.2500(10) 


-2.6681(10) 


-6.7476(10) 


1.9 


16.6259(18) 


7.9028(18) 


7.4211(18) 


1.2906(18) 




^MF 


unsmear 


.(i)MF 
APE 


[Mr) 

HYPl 


HYP2 




0.1 


-6.8662(15) 


-1.1655(15) 


-0.7654(15) 


2.1322(15) 




0.2 


-6.5422(15) 


-0.8342(15) 


-0.4337(15) 


2.4690(15) 




0.3 


-6.2962(15) 


-0.5805(15) 


-0.1796(15) 


2.7284(15) 




0.4 


-6.0971(15) 


-0.3732(15) 


0.0281(15) 


2.9419(15) 




0.5 


-5.9290(15) 


-0.1964(15) 


0.2054(15) 


3.1252(15) 




0.6 


-5.7847(15) 


-0.0430(15) 


0.3593(15) 


3.2855(15) 




0.7 


-5.6588(15) 


0.0926(15) 


0.4955(15) 


3.4285(15) 




0.8 


-5.5474(15) 


0.2143(15) 


0.6178(15) 


3.5580(15) 




0.9 


-5.4487(15) 


0.3240(15) 


0.7282(15) 


3.6761(15) 




1.0 


-5.3608(15) 


0.4236(15) 


0.8285(15) 


3.7848(15) 




1.1 


-5.2817(15) 


0.5151(15) 


0.9209(15) 


3.8860(15) 




1.2 


-5.2117(15) 


0.5985(15) 


1.0052(15) 


3.9799(15) 




1.3 


-5.1490(15) 


0.6756(15) 


1.0833(15) 


4.0683(15) 




1.4 


-5.0930(15) 


0.7471(15) 


1.1559(15) 


4.1522(15) 




1.5 


-5.0428(15) 


0.8141(15) 


1.2242(15) 


4.2327(15) 




1.6 


-4.9973(15) 


0.8778(15) 


1.2893(15) 


4.3112(15) 




1.7 


-4.9542(15) 


0.9408(15) 


1.3540(15) 


4.3908(15) 




1.8 


-4.9080(15) 


1.0089(15) 


1.4242(15) 


4.4774(15) 




1.9 


-4.8418(15) 


1.0994(15) 


1.5170(15) 


4.5888(15) 



Table 34. Numerical values of ^""^"^ and z^Vm'^ (DBW2). 
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, (pa) MF 



unsmear 



APE 



HYPl 



HYP2 



0.1 
0.2 
0.3 
0.4 
0.5 
0.6 
0.7 
0.8 
0.9 
1.0 
1.1 
1.2 
1.3 
1.4 
1.5 
1.6 
1.7 
1.8 
1.9 



8.8571(13) 
9.4715(11) 
9.76072(77) 
9.92818(67) 
10.03572(59) 
10.10844(54) 
10.15940(61) 
10.19583(64) 
10.22264(55) 
10.24404(58) 
10.26233(60) 
10.28122(78) 
10.30293(83) 
10.33363(89) 
10.37505(80) 
10.44470(87) 
10.56406(95) 
10.7954(13) 
11.3407(18) 



13.0847(13) 
13.6799(11) 
13.95403(78) 
14.10930(67) 
14.20697(60) 
14.27174(55) 
14.31641(62) 
14.34799(64) 
14.37123(55) 
14.39024(58) 
14.40716(61) 
14.42563(79) 
14.44776(83) 
14.47956(90) 
14.52260(81) 
14.59411(87) 
14.71514(95) 
14.9472(13) 
15.4907(18) 



13.5744(13) 
14.1649(11) 
14.43542(78) 
14.58812(68) 
14.68399(60) 
14.74759(55) 
14.79161(62) 
14.82300(65) 
14.84646(56) 
14.86603(59) 
14.88385(62) 
14.90351(79) 
14.92708(84) 
14.96057(90) 
15.00548(81) 
15.07897(88) 
15.20198(96) 
15.4358(13) 
15.9803(18) 



16.9632(13) 
17.5339(12) 
17.78944(80) 
17.93050(69) 
18.01739(63) 
18.07422(58) 
18.11334(64) 
18.14144(67) 
18.16305(58) 
18.18206(62) 
18.20048(64) 
18.22178(81) 
18.24793(86) 
18.28479(92) 
18.33365(83) 
18.41142(90) 
18.53858(98) 
18.7756(13) 
19.3207(18) 









(MMF) 






unsmear 


APE 


HYPl 


HYP2 


0.1 


55.6769(34) 


51.6386(34) 


52.5561(34) 


51.5545(34) 


0.2 


35.8730(16) 


36.7647(16) 


37.5663(16) 


39.4101(16) 


0.3 


29.0360(12) 


31.6064(12) 


32.3708(12) 


35.1870(12) 


0.4 


25.44278(95) 


28.88360(95) 


29.63050(95) 


32.95391(95) 


0.5 


23.1386(10) 


27.1312(10) 


27.8687(10) 


31.5163(10) 


0.6 


21.46410(79) 


25.85382(79) 


26.58617(79) 


30.46951(79) 


0.7 


20.13234(67) 


24.83618(67) 


25.56582(67) 


29.63781(68) 


0.8 


18.99324(66) 


23.96537(66) 


24.69400(66) 


28.92919(66) 


0.9 


17.95227(69) 


23.16938(70) 


23.89827(70) 


28.28431(70) 


1.0 


16.94561(81) 


22.40052(81) 


23.13075(81) 


27.66493(82) 


1.1 


15.91551(80) 


21.61484(80) 


22.34745(80) 


27.03545(80) 


1.2 


14.80158(97) 


20.76652(97) 


21.50263(98) 


26.35916(98) 


1.3 


13.5259(11) 


19.7962(11) 


20.5371(11) 


25.5886(11) 


1.4 


11.9728(10) 


18.6156(10) 


19.3631(10) 


24.6537(10) 


1.5 


9.9467(15) 


17.0763(15) 


17.8331(15) 


23.4370(15) 


1.6 


7.0623(14) 


14.8844(14) 


15.6548(14) 


21.7052(14) 


1.7 


2.4328(18) 


11.3644(18) 


12.1571(18) 


18.9236(18) 


1.8 


-6.5893(21) 


4.4999(21) 


5.3361(21) 


13.4952(21) 


1.9 


-33.2517(36) 


-15.8057(36) 


-14.8422(36) 


-2.5812(36) 



Table 35. Numerical values of 4^"^^^ and z^'"'^)^^ (DBW2). 
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